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The spectrum of cosmological perturbations in the context of the one-bubble open inflation model 
is discussed, taking into account fluctuations of the metric. We find that, quite generically,thin wall 
single field models have no supercurvature modes. However, single field models with supercurvature 
modes do exist. In these models the density parameter Q becomes a random variable taking a 
range of values inside of each bubble. We also show that the model dependence of the continuous 
spectrum for both scalar and tensor-type perturbations is small as long as the kinetic energy density 
. of the background field does not dominate the total energy density. We conclude that the spectrum 

' of the density perturbation predicted in the single-field model of the one-bubble open inflation is 

^\ , rather robust. We also consider the spectrum of scalar and tensor perturbations in a model of the 

Hawking- Turok type, without a false vacuum. 
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I. INTRODUCTION 



J> ■ In the standard inflationary universe scenario, both the flatness and the homogeneity problems are simultaneously 
C \ solved by the same mechanism, i.e., the accelerated expansion of the cosmic length scale. Therefore, if we want to 
solve the homogeneity problem, the flatness problem is automatically solved. This would make it impossible to create 
an open universe. Recently, however, attention has focused on a scenario which solves this difficulty. It is called the 
one-bubble open inflation scenario. The basic idea was first proposed by Gott III There are two main classes 

of models that realize Gott's idea: the single-field model, which was developed in ||,||], and the two-field models 
, proposed in S. Here we shall basically focus on the single-field model, but most of the results presented in this paper 
can be applied to models with two fields. 

In the single-field model we assume a scalar field with a potential such as the one shown in Fig. 1. Initially, the 
field is trapped in the false vacuum, where the vacuum energy drives the exponential expansion of the universe. This 
exponential expansion solves the homogeneity problem. After a sufficiently long period of inflation, bubble nucleation 
occurs through quantum tunneling. This process is described by the 0(4)-symmetric bounce solution 0, a solution of 
the Euclideanised equation of motion that connects the initial and final configurations. These configurations represent 
the state before and after tunneling, respectively. In the lowest WKB order, the classical evolution after tunneling is 
determined by the analytic continuation of the bounce solution. After the analytic continuation, the 0(4)-symmetry 
changes into the 0(3, l)-symmetry, which is just the symmetry of an open Friedmann-Robertson- Walker universe. 
Thus, in the lowest order approximation, we obtain a model that explains the creation of an open and homogeneous 
universe. However, at this stage the universe is almost empty and the cosmic expansion is dominated by the curvature 
term in the Friedmann equation. Hence, as long as only ordinary matter fields are assumed, the universe would stay 
curvature dominated forever. Therefore, in order to solve the entropy problem, a second stage of inflation inside the 
nucleated bubble is needed. 

Quantum 
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de Sitter expansion reheating 

FIG. 1. A schematic picture of the inflaton potential for a single-field model. 
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FIG. 2. A schematic picture of the potential of the second infiaton field for a two-field model. 



Although this model is simple, we need to assume a rather unusual form of the infiaton potential (see however ||). 
To solve this problem of naturalness, models with two fields were proposed, where two different inflatons would drive 
the first and second stages of inflation. The first infiaton field, a, has a double well potential. When the cr-field is in 
the false vacuum, the energy density of the universe is dominated by the vacuum energy of the cr-field. At this stage, 
the second infiaton field, (j>, has the potential shown by the dashed line in Fig. 2. After a sufficiently long period of 
inflation the cr-field makes a phase transition to the true vacuum through bubble nucleation. It is assumed that due 
to the coupling to the a field, the potential of the 0-field changes into the one shown by the solid line in Fig. 2. Then, 
the slow roll of </> from the old minimum to the new one drives the second period of inflation in the nucleated bubble. 

In some cases, the instanton describing tunneling has an approximate zero mode which can shift back and forth the 
value where the slow roll field lands after tunneling. In those cases, the density parameter in the resulting universe 
becomes a random variable, and one has an ensemble of large regions inside of each bubble where the density parameter 
n (measured at a fixed temperature) takes a range of values. This can be understood as follows. 

The existence of an approximate zero mode of the instanton implies that, in the real time evolution of the bubble, 
there is a supercurvature mode (whose wavelength exceeds the curvature scale). Quantum excitations of this mode 
can change the average value of the slow roll field inside the bubble on scales larger than the Hubble radius. As a 
result, the amount of slow roll inflation in different regions will be different, leading to a range of values of Q [Q. 
In particular, even if the instanton leads to a non-inflating value of the field after tunneling, localized fluctuations 
can create inflating islands which locally resemble open universes. This scenario is called "quasi-open" inflation |^,^). 
Quasi-open inflation was originally discussed in the context of two-field models, but as we shall see, models with a 
single field can also have supercurvature modes with similar consequences. 

Finally, Hawking and Turok |To[ have recently proposed that it is possible to create an open universe from nothing in 
a model without a false vacuum. The instanton describing this process is singular, and therefore its validity has been 
subject to question [jOj. Nevertheless, it has also been pointed out that the quantization of linearized perturbations 
in the singular background is well posed jl^ , |l3|| . Therefore, provided that one can make sense of the instanton by 
appealing to an underlying theory where the singularity is resolved, it appears that the details of that theory need 
not be known in order to calculate the spectrum of cosmological perturbations. 



FIG. 3. The Coleman-De Luccia instanton for the case that the spacetime dimension is two, embedded in the three 
dimensional Euclidean space. 
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FIG. 4. The analytic continuation of the Coleman-De Luccia instanton to the Lorentzian region. The lower half is the 
same as that shown in Fig. 3 and represents the Euclidean region. The upper half represents the continuation to the Lorentzian 
region. 

In order to construct an acceptable cosmological model, we must also show that the expected spectrum of primordial 
fluctuations is compatible with observations. Since there is an initial de Sitter expansion phase, it is natural to assume 
that the state before tunneling is the so-called Bunch-Davies vacuum state. Then the question that we must answer 
is the following: "What quantum state results after 0(4)-symmetric bubble nucleation starting from the Euclidean 
vacuum state?" In the lowest order WKB approximation, the tunneling process is described by the Coleman-De Luccia 
bounce solution |2|. In order to extract information about the quantum fluctuations, it is necessary to develop the 
WKB approximation to the next order [jl4| . The resulting prescription to find the quantum state after tunneling can 
be summarized as follows. 

First of all, since we shall deal with the metric perturbations which contain gauge degrees of freedom, we need to 
find the reduced action that contains only the physical degrees of freedom. We denote one of these degrees of freedom 
by ip(x). With an appropriate choice of the variable <p(x), the action can be reduced to that of the massive scalar 
field with variable mass, 
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gd x. 



(1.1) 



From this reduced action, one can calculate expectation values following standard methods. First, we need to find a 
complete set of functions «k which obey the field equation, 



0. 



(1.2) 



in the Lorentzian region, and which are regular on the lower hemisphere shown in Fig. 4. This restriction is due to the 
choice of the Bunch-Davies vacuum state before tunneling. However, the same restriction arises in the no-boundary 
proposal for the wave-function of the Universe |lq ], if Fig. 4 is interpreted as describing the creation from nothing of 
a universe containing a bubble. The same restriction will therefore be used in the case of the Hawking- Turok model 
discussed in Appendix D. Here, it is important to note that "a complete set" means a set of all modes which can be 
Klein-Gordon normalized on a Cauchy surface S of spacetime, 



/ dY,^ 1 " {vkd„v k > - (<9„«k)wk'} = <>kk' 



(1.3) 



Once they are obtained, the quantum fluctuations of the field are described by the "vacuum state" , | "J) , that satisfies 
Sk|^) — for any k where the annihilation operator Sk is defined by the decomposition of the perturbation field: 
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(w k a k + «ka k ). 



(1.4) 



Thus the mode functions play the role of positive frequency functions. 

The most tedious part in the above prescription is to find the reduced action. This is due to the fact that the 
time-constant hypersurface that reflects the maximal symmetry of background solution is not a Cauchy surface, 
and hence it is not appropriate to normalize modes. Therefore, we need to perform the reduction of the action on 
spacelike hypersurfaces which cut right through the bubble, and which are therefore not homogeneous. This renders 
the decomposition into scalar, vector and tensor modes into a rather unfamiliar form. In the end, however, the three 
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standard physical degrees of freedom are identified |16[ . When they are analytically continued to the open universe, 
one becomes the usual scalar-type perturbation and the other two becomes the even and odd parity tensor-type 
perturbations. 

Along the lines of the above prescription, the prediction of the power spectrum of perturbations was given by many 
authors |17] , [L8| , |2C|] . Until recently, however, the reduced action including the metric perturbations was not known, and 
for the scalar-type perturbation the studies were limited to the case in which the metric perturbation was neglected. 
Clearly, this should be a reasonable approximation in the weak gravity limit. Nevertheless, it turns out that the 
correspondence between the case where metric perturbations are neglected and the case where they are included is 
rather non-trivial. 

In jl6| we found that, at least formally, the metric perturbation has a dramatic effect on the spectrum for the scalar 
modes, even in the limit of weak gravitational coupling. When metric perturbations are ignored, there is a discrete 
wall fluctuation mode which gives a finite contribution to observables. However, when we incorporate the metric 
perturbation, this mode disappears. Fortunately, the contribution from the disappeared scalar-type perturbation 
reappears in the low frequency spectrum of even-parity tensor modes 20, 21[. In particular, we found that in the 
weak gravity limit the change of assignment of the wall fluctuation modes from scalar to tensor does not affect the 
observational quantities, such as the multipole moments of the temperature fluctuation of the cosmic microwave 
background. Note that the wall fluctuation mode is the Goldstone mode associated with the breakdown of spatial 
symmetry, and we find that it is "eaten up" by the gravitational degrees of freedom. This is reminiscent of the fate 
of Goldstone bosons in gauge theories. 

In the present paper we shall continue the work of |l6| ], investigating the full power spectrum of scalar perturbations. 
In particular, we shall consider the question of the existence of supercurvature scalar modes ]l8|,p^|. In the previous 
studies in which the metric perturbation was neglected it was found that, quite generically, there are no supercurvature 
modes for thin wall single-field models (other than the wall fluctuation mode). It was also shown that the contribution 
to the CMB anisotropy due to the continuous spectrum is almost independent of the details of the potential barrier. 
However, it is still uncertain whether such a robust prediction of the spectrum stays correct or not when we incorporate 
the effect of the metric perturbation. Here we shall clarify this point. As for the tensor-type perturbations, the studies 
in which the perturbation of the scalar field was neglected was done in Fortunately, this was found not to 

alter the resulting spectrum at all. However, the effect of the evolution after nuclcation of a bubble has not been 
studied so far. 

In previous studies, a simplified model was used to describe the background geometry inside the bubble, assuming 
that it was a pure de Sitter space. In this paper, we extend our analysis also taking into account that both the 
background geometry and the scalar field configuration are non-trivial. 

The paper is organized as follows. In Section II we first describe the lowest WKB order picture of the bubble 
nucleation. After that, we briefly recall the resulting reduced action for the perturbation on this background obtained 
in our previous paper |l6| |. In Section III we consider the question of existence of scalar supercurvature modes. For 
this purpose, several new techniques are introduced. In particular we show that the spectra in the weak gravity 
limit and in the case without gravity can be related to each other using the formalism of supersymmetric quantum 
mechanics. We also introduce a method for reconstructing the inflaton potential barrier starting from a given bubble 
profile for the tunneling field. This is used to find examples which do posess supercurvature modes. In particular, we 
discuss an generalization of the flat space Fubini instanton [ p7| to de Sitter space. This example is useful in order to 
illustrate the role played by the scalar supercurvature modes. In Sections IV and V we consider the spectrum of the 
scalar-type and tensor-type perturbations, respectively. Section VI is devoted to the conclusions. 

A number of issues are discussed in the appendices. In Appendix A we perform the canonical quantization of 
cosmological perturbations using the Dirac formalism. We recover the results of our earlier work |l6[ |, where the 
quantization was performed using the Fadeev-Jackiw reduction method. In Appendix B we derive some bounds on 
the "superpotential" which appears in the Schrodinger fluctuation operator for scalar modes. In Appendix C, we 
consider a class of models for which the evolution of perturbations can be solved exactly. Finally, in Appendix D we 
find the spectrum of cosmological perturbations in a soluble model of the Hawking- Turok type. We use the notation, 
k = 8irG, and adopt the units c = % = 1. 



II. REDUCED ACTION FOR PERTURBATIONS 



We consider the system that consists of a minimally coupled single scalar field, $, coupled with the Einstein gravity. 
We denote the potential of the field <& by V(Q). In order to describe the 0(4)-symmetric bounce solution, we consider 
a Euclidean configuration, 

ds 2 = a(rj E ) 2 [drj% + d X 2 E + sin 2 XE(d9 2 + sin 2 9dtp 2 )] , 
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$ = </>( Ve ). (2.1) 
Then, the equation of motion for the background quantities becomes 

<j>" + 2H(f>' - a 2 dV = 0, (2.2) 

H 2 1 = | Q^' 2 - « J VW) , (2-3) 

W'-W 2 + l = -^' 2 , (2.4) 

where the prime denotes a derivative with respect to jje, Ti '■= a! /a and dV := dV (<j>) / d<j>. We note that the above 
three equations are not independent. One of them can be derived from the other two. The bounce solution is the 
solution that satisfies the above set of equations with the boundary condition that </>' — ► 0, a — ► and H — + ±1 at 
j]e —> Too- The topology of the solution is the 4-sphere. Hereafter, we refer to scalar field and the scale factor of the 
background bounce solution by <f> and a. The two points at which a = are the centers of the 0(4)-symmetry. For 
convenience, we take t\e = 00 to be on the true vacuum side (i.e., inside the 0(4)-symmetric bubble), and call it the 
pole and that on the false vacuum side the antipole. 

Sometimes it is convenient to use the coordinate te defined by dTE — a(r]E)dr]E- In terms of it, the background 
equations are written as 

4> + 3-,j>-dV = 0, (2.5) 



a 

a \ 1 k 1 1 



,2 



- V{4>) , (2.6) 



a) a 2 3 V 2 

+ i = "?^ (2-7) 



where the dot represents a derivative with respect to te- 




FIG. 5. Conformal diagram of a de Sitter-like space with a bubble. 



Once we know the bounce solution, the background geometry and the field configuration in the Lorcntzian region 
are obtained by the analytic continuation of the bounce solution through the 3-sphere intersecting both poles. The 
coordinates in the Lorentzian region are given by 

Ve = Vc = -VR -7}i = VL + ^i, 
XE = -tXc + 2 = ~ l XL = -iXR , 

a = ac — iciR — ia^ ■ (2.8) 
Accordingly, the analytic continuation of the coordinate te is given by 
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cIte = drc = —idTR = idr^ 



(2.9) 



The coordinates with the indices C, i? and L cover regions C, R and L, respectively, in Fig. |5|. The region i? 
corresponds to the inside of the lightcone emanating from the center of the nucleated bubble, where a — and 
Vc = Ve — > +00. The region L is the inside of the lightcone emanating from the antipodal point where a — and 
•qc — * — co. The region C is the remaining central region. In region C, the metric is given by 

ds 2 = a(?/ c ) 2 [dr^ - dxc + cosh 2 Xc(d6 2 + sin 2 9d(p 2 )] . (2.10) 

One sees that surfaces that respect the maximal symmetry, i.e., the r\c =const. hypersurfaces are no longer spacelike. 
Instead of r/c, the coordinate xc plays the role of the time coordinate there. Note that the \C =const. hypersurfaces 
are not homogeneous. 

In order to find the quantum state after the bubble nucleation, we need to find a set of normalized mode functions 
that satisfies the regularity condition explained in Introduction. As mentioned there and can be seen from Fig. [^, 
no Cauchy surface exists in regions R and L. Thus, in order to obtain a set of normalized mode functions, we need 
to find the reduced action in region C, starting from the second order variation of the action including the metric 
perturbation on this background. This was done in our previous paper Jl6[ . There, we used the Faddev-Jackiw method 
for constrained systems to obtain the reduced action. In Appendix A of this paper, as an alternative and perhaps more 
familiar method, we present a derivation of the reduced action a la Dirac. Perturbations on a spherical symmetric 
background can be decomposed into even and odd parity modes. The odd parity modes do not contribute to the CMB 
anisotropy if we choose the position of the observer as the center of the spherical symmetry. Hence we concentrate on 
even parity modes. As mentioned above, we must work in region C, where the background configuration is spatially 
inhomogeneous. However, in region R or L, the background solution has the symmetry of the FRW universe. There 
the standard cosmological perturbation theory tells us that even parity perturbations for the Einstein-scalar system 
can be decomposed into the scalar and tensor-type perturbations. Therefore, one expects that even parity modes can 
be decomposed into two sets of decoupled perturbations even in region C and they are naturally identified with the 
scalar and tensor-type perturbations, respectively, when they are analytically continued to region R or L. This was 
shown to be true in Jl6[ . Hence, following the conventional terminology used in the cosmological perturbation theory, 
we call the corresponding modes the scalar-type perturbation and the tensor-type perturbation, respectively, even in 
region C. 

In this section, we work in region C where the Cauchy surface exists. Note that we have t\e = i]c and te = to- 
For notational simplicity, we omit the index C from the coordinate variables throughout this section. 



A. Reduced action for scalar-type perturbation 



First we consider the scalar- type perturbation. We introduce a gauge invariant variable Qs for the scalar- type 
perturbation which corresponds to the curvature perturbation in the Newton gauge in the open universe. Explicitly 
the metric and the scalar field perturbation in the Newton gauge is given as 



ds 2 = 1 + 25. dr 



4« 



- 2 ' a 2 



Hx 2 + cosh 2 X df! 2 ), 



Qi 



4/ca</> 



In terms of Qs, the reduced action is obtained in Appendix A as 



?(2) . 

5 sca — 



dr 



d 



cosh 2 x^-QtrnKs-H-Qfr 



K s {t(l + 1) + (K S -3)}Q 



(2.11) 



(2.12) 



where I and m are the eigenvalues of spherical harmonic expansion, Qg" 1 is the (£, m)-component of the spherical 
harmonic expansion, Qs =■ ^2 Qs ^m(^)) an d Ks is a derivative operator given by 



K s 



,3 12 



did na 2 4> 2 



dr a<p 2 Ot 

The relation of Qs to the variable q introduced in |l(| is 



(2.13) 



G 



Qs = 4:Ka(j>q . 

Putting q = £^ m Y^ m , we can use this relation to rewrite the reduced action ( [2.12| ) into the form, 



5 sca = J2\f d *f dr l 
where the operator O is given by 



cosh x— ~ — O 



Em 



dx dx 



qimQ + l) + (<5 - 3) C osh 2 x} q im 



(2.14) 



(2.15) 



O 



rf2 i k a,' 2 
^ + 2 



(2.16) 



The reduced action (2.15) exactly coincides with the one derived in [16j. Notice that in terms of Qs the action 
contains the operator Ks, which is not Schrodinger-like, whereas in terms of q it contains the operator 0, which is. 
Although this is of no essential significance, the second form is more convenient when normalizing the modes and 
when discussing the general form of the primordial spectrum. 

Now let us construct a set of positive frequency functions corresponding to the variable q. It is decomposed as 



(2.17) 



where A/jf is a normalization constant, and the sum is understood as a sum over the discrete spectra and as an integral 

over the continuum one. In this expression, q p {i]) is a spatial eigenfunction of the operator O with the eigenvalue 
p 2 + 4. That is, q p {rf) satisfies 



Oq p = 



d 1 K A' 2 



q p = (p 2 + 4)qP. 



The equation that determines the time evolution is found to be model independent and becomes 

1 



cosh 2 x ®X 



d ^ 2 d 

cosh %— 



" dx cosh 2 



X 



F l = (p 2 + i)f pt - 



(2.18) 



(2.19) 



Aside from the normalization factor, the positive frequency function is determined by the Bunch-Davies regularity 
condition, namely, it must be regular at x — (Xe = 7r /2)- The solution is 



f p \x) 



\ T{iy + £ + ljTj-ip + e+l) p -e-i/2 



F(ip)T(— ip) coshx 



jp-l/2 



\i sinhx), 



(2.20) 



where is the associated Legendre function of the first kind, and we have fixed the normalization factor so that the 
analytic continuation of Y ptm (x, f2) := / p '(x)y^™(^) to region R or L satisfies 



J sinh 2 x dx dtt Y plm YP' e ' m ' = 8{p - p')5u>5 n 



for p 2 > 0, where x = XR or XL- Also, these analytically continued harmonics satisfy 



(2.21) 



(2.22) 



where ^ 3 ^A is the Laplacian on the unit s paceli ke hyperboloid. For p 2 < 0, the normalization factor in Eq. ( 2.2C ) is 
not suitable. In this case, the solutions to (2.19) will be denoted by f A,i (x), where A 2 = — p 2 . The positive frequency 
functions f A ' l {x) ar e taken as the usual Bunch-Davies ones, 



,a,/, ^_ / r(A + i + i)r(-A + i + i) . 1/2 

2cosh x A - 1/2 



f A < e (x) 



(i sinhx), 



(2.23) 



and we define y ' (x,Q) ■= f (x)Y£m(ty- Notice that with this choice of positive frequency functions, y A < em are 
already Klein-Gordon normalized in the (2+1) dimensional sense, i.e. 
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i cosh x 



df 



Ox 



-f A ' e - f 



A.f 



df 



A.f 



dflYf m Yf'm' — 5m>5, 



U'<->mm' 



(2.24) 



These mode functions y A ' tm also satisfy Eq. (2.22) in region R 



In order to quantize the perturbations, it is convenient to write the action (2.12) in the canonical form (1.1). Intro- 
ducing a new variable q := £ + 4 a~ 1 J\f^q p (r})Y vtm (x, fi) +E V4 - A 2 a"W| Q' A (r?)y A ^ m (x, O), the action for 
q turns out to be the one for a scalar field with ^-dependent mass. Thus the normalized positive frequency functions 
are determined by calculating the ordinary Klein-Gordon norm for q. Namely we require 



j(p 2 +4)cosh 2 x |A^| 2 |^F-r^| 



dr] q p qP' 



2p(p 2 + 4) sinh 7rp |W S | : 



for the continuous spectrum (p 2 > 0), and 



(4-A 2 Mi| 2 



dr\ q p q p ' = 5{p — p ) 



drjq' q> A ' = Saa'j 



(2.25) 



(2.26) 



for the discrete one (p 2 < 0). Notice that the Klein-Gordon normalization reduces to a Schrodinger-like normalization 
for q p (once a convenient choice of A^f has been made) 



The variable Qs can be also decomposed in the same way as q. From Eq. (2.14), the corresponding spatial 



eigenfunction Q p s {t) is related to q p (r) as Q p s — y / ~Ka<j>q p . The equation satisfied by Q p s is given by the eigenvalue 
equation for the operator Kg as 



For later convenience, we also write down the equation for the scalar field fluctuation in the Newton gauge: 



(2.27) 



(2.28) 



where d 2 V := d 2 V ld<j> 2 and tpN has been expanded in the same way as q. We note that from the Newton gauge 
condition ( A28) in Appendix A, we have 



Q s = AKa(j){p p N = 4k— (acj)q p ) 
dr 



Using this, Eq. (2.28) may be rewritten as 



..p . Q a -P P ' + 

The analytic continuation of the variables Qs, and q to region R or L is given by 
Qs =iQs = *Qs> = <Pn = <Pn> 1c = -iQR = i( lL, 



(2.29) 



(2.30) 



(2.31) 



where the indices C, L and R denote the variables in the respective regions. The evolution equation in an open universe 
inside the bubble is given by the analytic continuation of Eq. (2.18), or Eq. ( |2.30 ) supplemented by Eq. ( [2.29 ). We 
note that the form of the equations does not change after the analytic continuation except for the signature in front 
of d 2 V in Eq. ( pC| ). 



B. Reduced action for tensor-type perturbation 



In this subsection, we consider the tensor-type perturbation. We introduce a gauge invariant variable Qt that 
corresponds to the tensor-type perturbation in the transverse-traceless gauge in the open universe. Let us denote the 
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metric on a unit 2-sphere by dVl 2 = &ABda A da B , and the covariant derivative and the Laplacian with respect to &ab 
by ( 2 )y and ^A, respectively. The tensor- type perturbation is expressed as 



ds = dr 



dx + [or + 



cosh 2 xJ V cosh 2 x 

+4 (9 x +tanh X ) (2) V A ( (2) A^Qt d X da A 

9 / - ( (2) A)-2 

+4 cosli x 2 - K T + ^ h tanh x<9 x 

V 2 cosh x 

x ( 2 >V^ 2 >V B - i<W 2 )A 
=: (is 2 -) + a 2 hij[QT\dx q 'dx\ 



cosh x <i" 



( (2) A( (2) A + 2)) -l gTdCT A dCT i 



where i, j = x, 9, <P and 



if T = 



_9_ 3_3_J_ 

9t <9t a 2 



(2.32) 



(2.33) 



The action for Qt is obtained in Appendix A and it is given by 



5', 



(2) 

GW 



E 

Z.m 



k{£ -!)£(£ +!)(£ + 2) 



d\ 



dr 
2^3 



cosh 2 x Q^K T (K T - l) Q*» 
ox dx 



-QpK T (K T - 1) {£{£ + 1) + K T cosh 2 x } Ot™ 



(2.34) 



where Q l ™ is defined by Qt ='■ Qfr'Ytm, as before. 

Again, the operator Kj- is not of Schrodinger type. Therefore, it is convenient to change to the variable w introduced 
in |lq], which is related to Qt by 



aP_d_Qr 
k dr a 2 



(2.35) 



Using this relation, the reduced action ( |2.34 ) can be rewritten as 

c(2! 



E 



2 dw hn <9«; 
cosh x — ~ — *v 



9x % 
{£(1 +l) + (K+l) cosh 2 x} 



to 



tin 



where the operator K is given by 



K 



92 j_ K A> 2 



(2.36) 



(2.37) 



This coincides with the reduced action obtained in Rig] . 

Now we construct a set of positive frequency functions. First we use the variable w. As before, we decompose w as 



(2.38) 



Note that the time-dependent part, f pi (x)i is t ne same as in the scalar- type perturbation. In this expression, the 
spatial function w p (rf) satisfies the eigenvalue equation for the operator K with the eigenvalue p 2 : 



Kw p 



d2 j_ k a> 2 
d^ + 2^ . 



w p = p 2 w p . 



(2.39) 



It is important to note that the spectrum of p is continuous with p 2 > because the potential term §^' 2 is positive 
definite and vanishes at infinity. Then, defining a new variable w := \J ij^ \Y(j+i)(t^ T) a~ 1 Afp l w p (ri)Y pem (x, Q), 







the action reduces to the one for a scalar field, as in the case of scalar-type perturbation. Therefore, the normalization 
condition becomes 



2 up 2 



2p sinh irp \ Af 



T|2 r oo 



{£-!)£{£ +!){£ + 2) 



dr]iv p wP' = 8{p — p'). 



(2.40) 



Since what we are interested in is the metric perturbation, we now consider Qt- It can be also decomposed as 
Qt = J2J\TpeQT P {T)Y pem (x, fi), where the spatial function Qt satisfies the eigenvalue equation for the operator Kt- 



KtQt = -a 4 



d 2 + g " d 
dr 2 a dr 



4 = (p 2 + i)Q t 



Then the metric perturbation hij [Qt] is given by 



11 _ jl/T Qr.y(+)ptm 

{£-!)£(£ +!)(£ + 2) pi a 2 11 



where 



Yi 



(+)pt m = 1 {£-!)£(£ +!)(£ + 2) 
' 2 V 2p 2 (p 2 + l) 



: h lj [Y plr ' 



(2.41) 



(2.42) 



(2.43) 



is defined so that its analytic continuation to region R or L gives the normalized even parity 



The tensor Y i 
tensor harmonics on the unit 3-hyperboloid [|23j, i.e. 



d 3 x^r'l jj 'Y}+ )pim Yf/ m ' = 5(p-p')5 u >6. 



(2.44) 



where jij is the metric on the unit 3-hyperboloid and 7 = det 7^ . Therefore it is more relevant to introduce a new 
normalization constant for the tensor perturbation: 



f/T . 2 / 2p 2 (p 2 + l)~ T 
J p -" V (£- 1)£(£ +!)(£ + 2) JVpt: 

with which the metric perturbation is expressed as 

In terms of NJ , the normalization condition ( |2.40] ) is re-expressed as 

Kpsmhirp\Np\ 2 



dr]W p wP' 



5{p-p). 



2ir(p 2 + 1) 

The relations between the analytically continued variables are given by 

Qt = Qt — Qt, w c = iwr = -iwl- 



(2.45) 



(2.46) 



(2.47) 



(2.48) 



III. SCALAR SUPERCURVATURE MODES 



The possible presence of scalar supercurvature modes (i.e. modes withp 2 < 0) has important physical consequences. 
If a model has supercurvature scalar modes then the density parameter Q becomes a random variable which takes 
different values in different places of the universe, and as we shall see, it may have a considerable spread if p 2 is close 
to —1 Also, if the spectrum of perturbations around a particular bubble contains a "supercritical" eigenvalue 

with p 2 < — 1 , this may indicate that we are not looking at the dominant channel for decay, and that there is another 
instanton with a lower Euclidean action (This is known as the "no supercritical supercurvature" mode conjecture 
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which has been shown to hold in several cases [E4| ) . In this section we shall analyze the conditions for the existence of 
scalar supercurvature modes in one field models of open inflation. We conclude that generically, these will be absent in 
the thin wall case, although they can be present for thick wall models (and even for thin wall models with sufficiently 
complicated potentials.) Intuitively, the reason is that for thin wall models the microphysical mass scales are large 
compared with the inverse size of the bubble, and hence all modes of fluctuation are "hard". In contrast, two-field 
models generically have supercurvature modes [jl8| even in the thin wall case. These are not due to the tunneling field 
but to the "soft" modes of the slow roll inflaton field - roughly speaking, the value of the slow roll field at the time 
of nucleation may be seen as an approximate zero mode of the action. 

First of all, in subsection [II A we shall investigate the weak gravity limit, comparing the scalar spectrum in the 
case when metric perturbations are included and the case when they are neglected. The formalism of supersymmetric 
quantum mechanics turns out to be very useful in this respect. Then, i n subs ection [II B we shall give a sufficient 
condition for the existence of scalar supercurvature modes. In subsection [II C, we describe a method to construct a 



model which has supe rcurv ature modes when the bubble wall is thick but does not have them when the bubble wall 
is thin. In subsection III D we consider the generalization of the Fubini instanton to a de Sitter background and use 
it to illustrate the role played by supercurvature modes in one field models. 



A. Supersymmetric Quantum Mechanics and the weak gravity limit 



For scalar perturbations the spectrum is determined by the Schrodinger operator O in Eq. (2.1S). This can be 
rewritten in the form 



OqP = [Co + («0'V2)]q p = (P 2 + 4)g p , 



where 



d 2 



(3.1) 



and we have defined 

W = <t>"/<t>'. 

From the boundary conditions satisfied by the instanton, it is easy to show that 

0' -> e ±2r > (ry -> Too) 

Hence, we have 

W^±2 (r)-> Too), 



(3.2) 



(3.3) 



and the potential term in (3.1) tends to 4 at infinity. Therefore the spectrum is continuous for p 2 > 0. Also, since the 
potential is not bounded below by its asymptotic value, there may be a discrete spectrum of supercurvature modes for 
p 2 < 0. The corresponding harmonics y A > lm have a wavelength longer than the curvature scale in the open universe 
(and are not normalizable on the t = const, hyperboloids [p5| .) 

Let us now consider the spectrum of perturbations of the scalar field $ when the metric perturbations are switched 
off. This is the limit where kinetic terms of the form K(f>' 2 are neglected, while the effective cosmological constant kV 
is kept finite. In that case, the background geometry is de Sitter and the equation for scalar field perturbations is 

V^Vju — d 2 V(j]c) ipk(i],Tc) = 0. Using the coordinates ( 2.1C ) and substituting the ansatz 
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— J n 



where Y p£m is understood as y A ^ irn for p 2 < 0, one readily finds the eigenvalue equation for % p : 

d 2 



th] 



2 X" 



[a 2 d 2 V + H' + H 2 + 3} X P = (p 2 + 4)% p . 



(3.4) 



(3.5) 



Using the scale factor for de Sitter space a(rf) — 1/ (H cosh 77) , this can be rewritten in the form 
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2 

drj 



2 X" 



[a 2 (d 2 V-2H 2 ) + 4)x p = (p 2 +4)x P 



Note that in de Sitter space the background field equation becomes 

1 



2 tanh 7] < 



H 2 cosh rj 



Taking the rj— derivative after multiplying it by cosh rj, we have 

' d 2 ( d 2 V~2H 2 

Izt-l- 

drj 2 



H 2 cosh rj 



dV(4>) = 0. 



<t>' = 0. 



It is clear from Eq. (3 
eigenvalue p 2 



that Eq. (3.6) has as its lowest eigenfunction the wall flu ctuation mode x 



.p=2i 



4. Hence, we can rewrite the operator in the left hand side of (3.6) as 



-ft 92 92 rxr/ 

V = h - — = 1- w 

drj 2 (/)' drf 



W z 



(3.6) 



(3.7) 



(3.8) 



with 



(3.9) 



where W is defined as in (3.1). In going from (|3J]) to (p.9|) we have used that the background metric is de Sitter. 



However, for future reference, we take (3.£) as the definition of the operator V , even in the case of strong gravity. 

Thanks to the formalism of supersymmetric quantum mechanics p6| , the correspondence between the spectrum 
of V and that of O in the limit of weak gravity becomes rather straightforward. It turns out that Oo and V are 
supersymmetric partners of each other, so neglecting the term k^' 2 , both spectra coincide (except for the ground 
state). Indeed, in terms of the "supersymmetry" generators 



Q = d v + W, Q f = -dr, + W, 



we have 



Oo = Q^Q, V = 



Note that if supercurvature modes are absent for the operator Oo, they are absent for O since the neglected term 
K<f>' 2 /2 is positive definite. Normalized eigenstates x V s °f T 3 (which we may call "bosons") can be transformed into 
normalized eigenstates q p s of Oq ( "fermions" ) and vice versa by acting upon them with the SUSY generators: [] 



QV s - (P 2 



4)^5, 



(3.10) 



Qq p s = (p 2 + 4) 1/2 X P S - (3-11) 
Note that, for convenience, in the above equations we have taken the eigenstates to be normalized as in the corre- 



sponding Schrodinger problem [and not as in (2.25)], i.e. 



9s«s=W» / XsXs= s p,p' ( 3 - 12 ) 

Therefore, the spectra coincide except for the ground state of V with p 2 + 4 = 0, whose eigenfunction Xg -2 * is 
annihilated by and has no fermionic partner. The reason is that this partner would satisfy the equation Qq s = 0, 



which due to the boundary condition (3.3) has no normalizable solutions. This is why the wall fluctuation mode 
disappears from the scalar spectrum when gravity, even if weak, is turned on. 

Where does the wall fluctuation mode go? Even in the absence of gravity, we may choose to eliminate the mode 
X p=2t from the scalar spectrum. Through a coordinate transformation, it can be rewritten as an even parity tensor 
mode with p 2 = plf , which is precisely the bottom of the continuum spectrum of the operator K corresponding 



1 The subindices S are included to indicate that these are eigenstates of the SUSY related operators Oq and V, and not of 
physical operators such as O. 
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to tensor modes, denned in Eq. (2.39). When the gravitational coupling is switched on, it turns out that the effect 
of wall fluctuations, which was concentrated in the p 2 = tensor mode, spreads to the long wavelength region of the 
spectrum, with eigenvalues M 



P 



2 < 



K(j)' 2 di] < 1. 



ft has been shown that these long wavelength tensor modes reproduce the effect of the wall fluctuation mode on the 
cosmic microwave background in the limit of small K(f> 12 pCj ]. 

In |lq| it was shown that, neglecting metric fluctuations and under rather generic conditions for the bubble profile 
in the thin wall limit, the operator V has the wall fluctuation mode as its only supercurvature mode. Specifically, 
the discussion given there assumes a step function model of d 2 V that has the following behavior of d 2 V . First, 
d 2 V — 2H 2 > on the false vacuum side followed by a thin region of d 2 V — 2H 2 < 0. Then there is a thin region of 
d 2 V — 2H 2 > followed by a region of d 2 V <C 2H 2 on the true vacuum side.[] It follows from our previous discussion 
that under those conditions the operator O will not have any supercurvature modes. 

We can also argue another case for the non-existence of a supercurvature mode as follows. Let us consider Eq. ( |3.6| ) 
for x p - We assume the case when the background solution is such that a region of /7 e //(ry) :— (d 2 V — 2H 2 )a 2 < 
exists only for a finite interval of 77. Let r\i (i — 1,2) be the two zero points of C/ e //! U e ff(r]i) = C^ e //fe) = 0. We 
assume 771 < 772- Now consider a one-node solution of x p ln t ne interval rji < rj < 772 with the boundary conditions 



X p/ ( T 7i) = X P '(V2) = 0. Let p 2 be the eigenvalue of this solution, 
a one-node solution x Po with p\ < for the original problem, x p 

771 < 77 < 772- Then it follows th at p 2 > p 2 . In other words, there is no supercurvature mode if p 2 > 0. 



Since U e ff is positive elsewhere, if there exists 
must be more warped than x p * m the interval 



As we shall see in subsection [II C , there exists a method to reconstruct a potential V from a given <f>' as a function 
of 77 provided it satisfies the boundary condition ( |3.2j ) and it is positive (or negative) definite. Let ^'(77) = 4>'(ar]) 
where a is a constant within the interval 771 < 77 < 772 and a — > 1 for 77 — > ±00. Then one can always construct a 
potential V from <fi' . Now let us consider a constant scale transformation 77 — * arj to the equation for x p within the 
interval 771 < 77 < 772: 



(p 2 +A)x p {ar,) = 



if 



1 



d 2 



d{ari)' 2 



4>'{ar]) d(ai]) 2 ^ ^ ' 



X p (uri) 



P 2 



drj 2 

[ -x p (v), 



1 d 2 

4>'{ri) d i 



x p (v) 



(3.13) 



where X P ( 7 ?) = X p ( ar l)- Thus we see x P is a solution if x p is so, and the eigenvalue p is given by 

p 2 +4 = a 2 {p 2 +4). 



(3.14) 



Then, taking x p to be the one-node solution with the eigenvalue p 2 constructed in the above, the corresponding 
one-node solution for the potential V has the eigenvalue 

Pl=a 2 (pl+4)-4. 

For a sufficiently large a, we have p 2 > and hence there will be no supercurvature mode for V. Since a 3> 1 
corresponds to taking the thin wall limit, this implies the supercurvature mode disappears in the thin wall limit for 
the class of potentials which are related by the above scale transformation. 

Unfortunately, neither the arguments given in |l^] , nor the above scaling arguments are completely general, and we 
have not been able to give a proof that thin wall always implies that there are no supercurvature modes. Therefore, 
in principle, one has to analyze each case separately. Nevertheless, the results given in the following subsections may 
provide some valuable guidance in this respect. 



2 It can be easily seen that letting the region of d 2 V — 2H 2 > thick or letting d 2 V on the true vacuum side larger will only 
make the existence of a supercurvature mode more difficult. 
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B. A sufficient condition for the existence of supercurvature modes and an SWKB estimate 



An upper bound on the lowest eigenvalue of the scalar operator O can be obtained by the "variational" method. 
Writing 

d = v -2W + 

we can readily evaluate the expectation value of O in the "test" quantum state whose wavefunction is given by <j)' 
(i.e. the true ground state of the operator V). This expectation value must be larger than or equal to the lowest 
eigenvalue p of O, which implies 



Therefore, a sufficient condition for the existence of a supercurvature mode (p 2 , < 0) is that the bubble profile satisfies 

dr) < 0. (3.15) 



A" 2 J/2 



Typically, the third term in the integrand can be neglected because it is Planck scale suppressed. Since the first 
term carries two more derivatives than the second, this indicates that we may expect supercurvature modes when the 
bubble walls are thick. 

Aside from this sufficient condition, we can also estimate the number of bound states of O using the supersymmetric 
version of the WKB quantization condition (SWKB) |2{|: 



/ \Jp 2 + 4 - W 2 dr\ = irn, 

Jb 



(3.16) 



where n = 0,1,2,... for the spectrum of the operator V, and n — 1,2,... for the spectrum of the operator Oq. 
The integration runs over the range where W 2 < p 2 + 4.^| Hence, the approximate condition for the existence of a 
supercurvature mode of Oq is that 



/ = / ^4 - W 2 di] > 7T 

Jb 



(3.17) 



As shown in Appendix B, this condition is not satisfied when the bubble walls are thin, so we do not expect super- 
curvature modes in this case. 



Unlike the sufficient condition (3.15), however, the condition (3.17) is the result of an approximation, and so it 



cannot be used as a rigorous criterion for the existence of supercurvature modes. 

C. A simple method to construct a model 

In this subsection we point out that for any given profile <p{vi) which is monotonous and has appropriate boundary 
conditions at infinity, we can always "reconstruct" a potential V{<p) for which (j) is a bounce solution. In particular, 
this freedom can be used to construct models which do have supercurvature modes. 



infinity. Then, the solution of (2.4) at rj — * ±oo given by 

n-i 



Let us consider that </>' is given. It can be chosen arbitrarily, provided that it has no nodes and satisfies (3.2) at 

= C±e 2r >. 



H + l 



" K '" ' * 1 " "-• ''■ * 4 = 0. In this case the "turning points" are coincident and the 



Note that for the ground state n = we know that p 
SWKB approximation gives the exact answer. As noted in 



2q] , this approximation gives accurate answers at both ends of the 



allowed range of n (n = and n — » oo), so one may expect better results than in the usual WKB approach. 
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This solution, which has been obtained from taking <j>' strictly to zero, automatically satisfies the required boundary 
condition, TL — > =Fl as rj — > ±00. Once we know 7i, we can determine the scale factor by integrating Ti. Here an 
integration constant appear s wh ich will determine the overall scale of the potential (say the potential at the false 
vacuum Vf)- By using Eq. (2.6), the potential is determined as a function of rj 



V 



6-6H 2 - k. 
2kcl 2 



s' 2 



(3.18) 



Since by assumption <p(rj) is monotonic, we can express rj in terms of <fr and the shape of the potential V((f>) is uniquely 
determined up to the overall factor mentioned above. 

Here, we note that the positivity of the potential determined in this way is not guaranteed in principle. However, 
integrating the equation 



V 



(3.19) 



and using the boundary condition 4> — > as 77 — > ±00 , we have AV = 3 a 2 TL4>' 2 dri, where AV is the difference 

between the heights of potential in the false vacuum (rj — > —00) and at the nucleation point (77 — > 00). If K<fi' 2 is 
sufficiently small, then AV/Vp becomes small, and hence V stays positive definite. 

Having shown that V can be constructed for given <f>', let us now turn to the eigenvalue equation for the scalar- type 
perturbations in terms of </>, Eq. (2.18). If there is an eigenfunction with negative p 2 , the solution with p 2 = must 



have, at least, one node. Hence, it is sufficient to examine the equation for q := q p °, which is written as 



drj 2 



k a> 2 
2^ 



U(r,) 



where 



U{rj) = -W + W 2 



(3.20) 



(3.21) 



Since the change in the overall amplitude of <f>' does not alter the function W, we can construct a model in which the 
second term is negligibly small without changing the potential U. Hence we neglect this positive definite term. As 
mentioned above, the positivity of V is guaranteed in this limiting case. 

Instead of specifying the mod el b y giving 0, we can specify it by giving W , which can be chosen arbitrarily as 
long as the boundary condition ( |3.3| ) is satisfied. Then, by choosing W to stay near zero for a sufficiently long time 
interval, U will have a wide negative- valued region and the solution for q = q p=a will have a node. On the other hand, 
if the wall is sufficiently thin, W' will dominate the potential U when \W\ is small. Then U may become positive 
definite, and in this case q would have no node. 

For instance, let us consider that W is given by 



W = -2tanh((?7-£)/A), 
where the parameter £ indicates the position of the bubble wall. Then U becomes 



1 



cosli (??-£)/ A 



(3.22) 



(3.23) 



When A > 1/2, the potential U becomes negative. Then q has a node, and hence there are supercurvature modes. 
On the other hand, when A < 1/2, the potential U is p ositiv e and there are no supercurvature modes. Of course 
these results are consistent with the sufficient condition (3.15), which for the present case can be shown to imply 
A > 3/4. Also, it is known that for the "superpotential" (3.22) the SWKB condition is exact, so (3.17) is equivalent 
to the condition A > 1/2. Note that A ~ 1 separates the thin and thick wall regimes, since for A <C 1 we may write 

wall thickness ~ a(£)A ^ a (0 ~ W& 11 radius. 



4 Note that the equations of motion for the background have the following symmetry: a global factor in the potential can be 
reabsorbed into the scale factor, and so its only effect is to change all physical distances by a constant factor. 
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The form of the potential V(4>) corresponding to (3.19) cannot be given in closed form for arbitrary A. However, as 
noted in for the case A = 1 it corresponds to the inverted quartic double well potential, which is discussed in 
the next subsection. 

Here we note that the above method for constructing a model does not prove that (j>(rj) will actually represent a 
relevant instanton, since for the same potential V there may be another solution which gives a larger decay rate. This 
issue is discussed in Q . 



D. Generalized Fubini instanton and the implications of supercurvature modes 



In flat space, the inverted quartic potential 

V f = -(A/4),/, 4 (3.24) 



has the following solution, known as the Fubini instanton |27 , 28 



^ at = VA^T7 (3 ' 25) 

Here, r is the distance to the origin of polar coordinates, and the arbitrary parameter p can be interpreted as the size 
of the instanton. This is therefore a one-parameter family of instanton solutions which interpolate between the "top" 
of the potential <fi = at r — > oo (which can be considered as "false vacuum") and an arbitrary <j) = (8/A) 1 / 2 /? -1 
down the hill at r — > (which can be thought as true vacuum). Note that in this flat space model there is no classical 
barrier between those two values of the field, and hence the instanton represents tunneling without barriers pSj]. The 



Euclidean action for ( 3.25 ) is 



Se = |£. (3.26) 
The fact that the solution exists for all si zes p and that moreover the action is independent of p is a consequence of 



the conformal invariance of the potential (3.24), which contains no mass scales. 

Since de Sitter space is conformally flat, we can generalize the Fubini instanton to de Sitter space. Therefore, for 
simplicity, we shall neglect the gravity of the bubble and we shall work directly on the 4-sphere. The metric can be 
written as a conformal factor times the metric of flat space 

s flat 

where 



ds 2 = A 2 ds 2 flat = A 2 [dr 2 + r 2 dS 3 ] , (3.27) 



H^ 1 is the radius of the 4-sphere and dS 3 is the metric on the 3-sphere. 
Under conformal transformations 

ds 2 flat ^ds 2 =A 2 ds 2 flatl <f>fiat —></> = A _1 ^// t (3.29) 

the Klein-Gordon operator has the following property (see e.g. p9|) 

[a-2H 2 }^ = A- 3 a flat flat . 

It follows that 

is a solution of the field equation □</> = dV/dcj) for the new potential 

V = 2H 2 <? - + const. = -^ 2 - v 2 ) 2 , (3.31) 
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where 



2H 2 



Therefore, the solution ( 3.30 ) generalizes the Fubini instanton to de Sitter space. 

Note that we have added a conformal mass term to the potential V{4>). The effect of this mass term is that now 
the potential looks like an inverted double well, with a local minimum at <j> = 0. The solution (|3.30[ ) now interpolates 
between 



vp 



(3.32) 



at the "pole" (or nucleation point) r = and 



vp 



at the "antipole" (or antipodal point) r — > oo. Here we have introduced the parameter p := pH/2. The generalized 
instanton no longer represents tunneling without barriers. It interpolates between a point inside the local well near 
4> = and a point outside this local well, beyond the local maxima at <j> = v. Therefore this is a Coleman-De Luccia 
instanton of the usual sort. For p = 1, both turning points coincide and we have <Pn = 4>a — v, so the Coleman- 
De Luccia instanton degenerates into a Hawking-Moss instanton, whe re th e field is constant t hroug hout the sphere. 
Introducing the variable r) = — \n(rH/2) the metric takes the form (2.1) and the instanton (3.30) takes the form 



cosh 77 
cosh(ry — £) 



(3.33) 



Here we have introduced a positive parameter £ which can be interpreted as the position of the bubble wall in the 
conformal coordinate r). It is related to p through 



The solution in the form (3.33) was given previously in p4| , where it was found using the method described in 



subsection IIIC Indeed, it corresponds to the case A = 1, which as mentioned there has a supercurvature mode. 

Let us now discuss the dynamics of the resulting bubble. Because of conformal invariance, the value of the Euclidean 
action in curved space is also independent of the value of p. Naively, one would then expect different types of bubbles 
to nucleate with equal probability, characterized by a different initial value from which the field starts rolling down. 
However, as we shall see, this picture is too simplistic. What actually happens is that fluctuations of the supercurvature 
mode cause different regions inside the same bubble to start from a different value of 4>n, so that all possibilities are 
sampled within just one bubble. 

Let us recall the form of the supercurvature mode in this model. Since ( 3.33| ) is a solution for any value of the 
parameter £, the associated zero mode 



(3.34) 



must be a solution of the linearized equations of motion about the solution (3.33). This perturbation onl y de pends 
on ?7, and not on the spatial coordinates on the hyperboloid. Hence, it must correspond t o a s olution of (3J3) with 



(E-22D 1 



-1 = -A 2 , 
Indeed, 



0, since only in this case the harmonic 3^ A ' 00 in the expansion (3.4) is constant [see Eq 



oc a(j])tpo oc 



— sinh(?7 — £) 
cosh 2 (ry — £) 



(3.35) 



solves Eq. (3.6) and it is, moreover, normalizable in the Schrodinger sense. However, this mode is real and hence 
it is not normalizable in the Klein-Gordon sense (its Klein-Gordon norm is zero). The reason is that its expected 
amplitude is actually infinite. Physically, this corresponds to the fact that all values of £ (that is, all values of the 
field <f>N at the time of nucleation) are equally probable. To handle this case, the zero mode has to be quantized in 
terms of "position" and "momentum" operators associated to the collective coordinate £, and not in terms of creation 
and annihilation operators accompanied by the corresponding Klein-Gordon normalized modes. For our purposes, 
however, it will be more instructive to break the degeneracy by perturbing the potential (3.31), so that not all values 
of <pN are equally probable. 



17 



Following pj], we add the term 



AV = eV 1 (4>), 



(3.36) 



to the potential (3.31), where e is a small parameter. In this case, the Euclidean action evaluated on the field 
configuration cf> = (f>^ + A(f>, where A(f> is a small perturbation of 0(e) which is otherwise arbitrary, is given by 

S(0 := AS E [<t>e + M = S%[<j> 6 + A0] + 2ir 2 e f V^(^) + 0(e 2 ). 



(3.37) 

Since cf>£ is a solution of the unperturbed potential, the field perturbation A(f> does not appe ar in the right hand side 
of this equation to linear order. Thus, to the lowest ord er in e, the instanton is still given by (3.30), but no w wit h the 
value of £ = £° that minimizes the right hand side of ( 3.37 ). This selects a "preferred" value ^ . Using ( 3.33|) , one 
can change the variable of integration in (3.37) from r\ to <j>£. We have 



Wir 2 



Q>-$)0.-p$)$- i v 1 (<i>)d4 ) 



where 
have 



ff 4 (1-p 2 ) 3 

ly is the dimensionless field. For definiteness, let us take V\ — if 4 4 (ao + a\4> + a2<fi 2 )- In this case we 



S(p) = Air 2 e 



ai 1 



Ct!2 1 



^41 



const. 



This expression has extrema at 



p = p = P±W 2 -l) 



1/2 



(3.38) 



where (3 = —5ai/12a.2. Of course, ( [3.38 ) is only valid for (3 > 1. F or < 1 the extremum is at p = 1 and the 
Hawking-Moss transition is the dominant channel. The double sign in ( 3.38 ) corresponds to the two different turning 
points for a given instanton, so actually there is just one extremum for each value of (3. 
One can estimate the spread in 0jv by using the formula for the probability of nucleation 



>N 



(3.39) 



where the exponent is calculated by substituting the function (f>^ which corresponds to a given value 4>n in ( |3.37| ). We 
have 



Using 



f'N = 4>Nd£, the variance of 4>m is given by 

<(A<M 2 > = 



d 2 S 



(3.40) 



In |24j it was shown that by adding the perturbation (3.36), the eigenvalue at p 2 = — 1 shifts to a new eigenvalue 

3H 2 d 2 S 



7 := P + 1 = 



4ir 2 v 2 dS, 2 



Introducing this expression back in (3.40) we have 



<(A<M 2 > = 



7T 2 7 



(3.41) 



Here po = 2poH 1 is the size of the bubble. Note that the eigenvalue 7 is of order e, and so as expected the spread goes 
to infi nity as we switch off the perturbation. Although the result (3.41) has been obtained here rather heuristically 
from (||9|), it can be justified more rigorously by studying the r.m.s. fluctuation of <f> in the 0(3,1) invariant quantum 
state, as it was done in || for two field models. In particular, by calculating the two point function it can be 
shown that the spatial correlations in </> at the time of nucleation decay on the co-moving scale r ~ 7 -1 , where r ~ 1 
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corresponds to the curvature scale on the open sections. Therefore, we have the picture that inside of each bubble 
there is an ensemble of regions of size r ~ 7 -1 where the value of 4>n is coherent. But the different regions of the 
ens emble take all possible values of 4>n, with statistical average given by the peak value 4>° N and with r.m.s. given 
by ( 3.4C ). The situation here can be contrasted with the one for the two-field models considered in ||. There, the 
mean value of the slow roll field <fi selected by the instanton solution is actually zero, at the bottom of the slow roll 
potential, and inflation is only due to localized fluctuations with r.m.s. given by the quantum state. Here, the mean 
value of 4> is nonzero, at a certain value (jP N from which the field starts rolling down. 

Also, one can show that r.m.s. amplitude of the supercurvature modes with I > is of order p^ , and not of order 
H as it is for subcurvature modes. This is a rather general conclusion, and physically it corresponds to the fact that, 
since these modes "live" in the bubble wall, they are excited due to the bubble acceleration, and not just due to 
the expansion of the universe. To avoid too large CMB anisotropics caused by the supercurvature modes (note that 
Pq 1 > if), in one field models we must 'tune" the parameters to some extent, so that the radius of the bubble is 
comparable to the Hubble radius. In our case, this is achieved for j3 ~ 1. 

Unfortunately, the inverted quartic potential is not too appropriate from the point of view of slow roll inflation. 
After nucleation the field would roll down to infinity in a finite proper time and the kinetic energy would soon dominate 
over the potential. Nevertheless, it appears that in any model with supercurvature modes the picture described above 
should be qualitatively similar. The different values of <j> & t the time of nucleation would then result in different 
lengths of the slow roll inflation, leading to an ensemble of regions with different values of the density parameter. 
This is the same picture that arises in the case of two field models studied in j?],^. More general single-field models 
with supercurvature modes may be viable phenomenologically, although they are likely to be rather fine tuned. This 
issue is currently under investigation. 



IV. SPECTRUM OF SCALAR-TYPE PERTURBATION 



A. The quantum state at t/r = —00 



In our previous paper |l6j , we showed that the general features of the spectrum of the tensor- type pertur batio n can 
be understood without specifying details of the potential model. This fact will be briefly reviewed in section V A . Here 
we show tha t an a nalogous argument also holds for the scalar type perturbation. For this purpose, it is convenient to 
use the Eq. (2. IS). We expand the quantum operator q corresponding to the variable q as 



Q = E 4t^V(^™( x , ^)+E a$ m M s A q> A (ri)y 



A,lr, 



( X ,0) + (h.c). 



(4.1) 



a (s) ~ (s) 

where a pem , a\\ m are the annihilation operators. Below, we shall analyse the behaviour of these modes in region C. 
The language of one dimensional scattering off a potential barrier will be used. Then we will analytically continue 
the modes inside region R and compute its evolution till the end of inflation. 

Noting the limiting behaviors 4>' cx e T2r,c at r\c — > ±00, the potential term of Eq. (2.18) tends to 



K 

r 



and Eq. ( 2.18 ) reduces to 



— < 



4 + 0(a 2 ) (nc -►±00), 



(Vc 



±00) 



(4.2) 



(4.3) 



Hence, for the continuous spectra one can choose the plane waves e Tipvc incident from rjc — ±00 as the two orthogonal 
solutions for q±. These incident plane waves interact with the potential and produce waves reflected back to r\c — > ±00 



with reflection amplitudes g±, and waves transmitted to r\c 
the limit r/c — > ±00 are given by 



=Foo with transmission amplitudes a±. That is, q± in 



and 



iql 



g+e 
<j + e 



ipvc ^_ e —ipvc _ 

ipVC 



<j- e 
Q- e 



iprjc 
-ipVC 



awnc 



(vc 



(Vc 



+00), 
-00), 



+00), 
—00). 



(4.4) 



(4.5) 
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Here we note that q + p and q_ p are not independent of q p + and q p _. Hence we restrict p to be non-negative. The 
reflection and transmission coefficients satisfy the relations 



|a + | 2 = l-|g + | 2 , |a_| 2 = l-|0_| 2 , 
<7_ =cr+, <J+Q- + = 0, 



(4.6) 
(4.7) 



which follow from the fact that the Wronskian of any two solutions of ( |4.3| ) is constant. It is easy to see that 
drj q p a q p a , = 2ir5(p — p')8 aa > . Thus, from Eq. ( 2.25 ), the normalization constant Afp is determined as 



A/" p s = (4p(p 2 + 4) sinh np) 1/2 . (4.8) 
The analytic continuation to region R of Eqs. (|4.4|) and (JO) gives the continuous modes in the limit tjr —>■ — oo as 



(4.9) 



q P R+ = e^ 2 Q + e- ipnR + e^^e*™*, 



1 P n- 



= 7rp/2 



with the constraint on the coefficients, \g+\ 2 + |c_| 2 — 1. These play the role of the positive frequency mode functions 
in region R. 

For the discrete part of the spectrum, the normalizable solutions q A can be choosen to have the form 



(rj -> too) . 



The normalization constant Af^ is determined from Eq. ( 2.2(i| ) , 

Ka= f(4-A 2 ) f" dv\q 

\ J — oo 



A|2 



-1/2 



The analytic continuation to region R of Eq. (4.1C) gives the discrete modes in the limit rjn — > — oo as 



A _ e iA7r/2 e A7) R _ 



(4.10) 



(4.11) 



(4.12) 



B. General formula for the power spectrum 



The important quantity that determines the primordial density perturbation spectrum as well as the large angle 
CMB anisotropics is the curvature perturbation on the comoving hypersurface, 1Z C . The comoving hypersurface is 
the one on which the scalar field fluctuation ip vanishes. In the present case, 1Z C is given by 



K p = -AC 



■s t *<t>'<i p 

2a 



H 



P 



Using Eq. (§]||, it can be written in terms of q alone, 



n p 



-AT" 



K(j/_ _H d_ 

2 a a<j)' 2 drjR 



(4.13) 



(4.14) 



The evolution of the scalar-type perturbation is governed by the analytic continuation of Eq. (2.18) to region R. It 
takes exactly the same form as the one in region C: 



dr) 2 R 



(4.15) 



Now, we will define a function that will gives us the value at the end of inflation of the curvature perturbation 1Z P . 
which is generated by a scalar per turba tion q p having the form e lpT,R on the lightcone when it enters inside region R 
(i.e. when r/R — > — oo). From Eq. (4.14), this "transfer function" for the scalar curvature perturbation is given by 



T p 



lim 

VR— ^end 



n 



k£ 

2 a a<fi' 2 drjn 



(4.16) 
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where, accordingly, q p is the solution of Eq. ( 4.1 5| ) which behaves as q p — » e lpr ' R when rjn — > — oo, and where we have 
drop the normalization constant A/If. Here ?7 cn d is the value of the conformal time tjr for which a — > oo. Taking into 
account that the quantum state for q p inside region R is given by Eqs. (|4.9|), it is easily seen that the continuous part 
of the primordial spectrum for 1Z C at the "end of inflation" can be computed from 



(l^l 2 ) = W + \ 2 + I^M 2 - 2|Af \ 2 \T P \ 2 (coshnp + Zi 



Using Eq. (4^8) we have 



(4.17) 



|T P | 2 / T p 

2^ + 4)sinh 7 rJ C ° Sh7rp + 5R ( g+ 7f 



(4.18) 



For the discrete modes we define the transfer function as follows, 



lim 



K<j>' H d 
2 a a(f>' 2 drjR 



(4.19) 



where q' A is the solution of Eq. (4.15) which tends to e Ar,R on the lightcone inside region R. The discrete part of the 
primordial spectrum for the curvature perturbation is then given by 



A|2\ 



,A|2 



Wl\ 2 \r s 



A,2 



(4.20) 



where is given by Eq. (4.11) 



Thus, in order to find the power spectrum for scalar- type perturbations one has to solve two separate problems. 
First, for the continuous spectra, one has to solve the scattering problem posed in the previous subsection in order 
to find the reflection coefficient p+ . Roughly speaking, this accounts for the scattering of scalar modes off the bubble 
wall, which determines the initial quantum state for these modes at the beginning of the second stage of inflation 
inside the bubble. For the discrete spectrum, one has to find the normalization of the modes, which will give us the 
"size" of their effect. Then, one has to compute the transfer function Tg, which accounts for the evolution of the 
modes during the second stage of inflation. 

The expression (4.18) is of general validity. Typically, however, both the reflection coefficient and transfer function 



can only be found numerically. In the appendices we discuss some exactly soluble models, but for the rest of this 
section we shall consider the weak backreaction case, which can also be treated analytically. 



C. Weak backreaction case 



We now consider the evolution of the scalar-type perturbation inside the bubble in the weak backreaction case. For 
notational simplicity, we set t — tr and omit the suffix R from all the variables in the rest of this subsection unless 
ambiguity arises. 

In general, one cannot assume the slow roll motion of (f> immediately after nucleation. Nevertheless, it can be assumed 
that the back reaction of the motion of <j) to the cosmic expansion is negligible at this stage when the universe is 
curvature-dominated. Then, after a few Hubble expansion times, the universe starts to expand exponentially and the 
slow roll assumption becomes valid. Hence it is fair to assume 

«tf 2 ^%^H 2 + \, (4.21) 

where H 2 := kV/3. We call this the small back reaction approximation, since it corresponds to the limit in which the 
back reaction of the inflaton field dynamics to the geometry is small. In this limit, a can be approximated by that of 
de Sitter space as long as the time duration considered, At, is not much greater than if -1 . Also, we shall make the 
approximation that d 2 V varies slowly compared with the expansion time. Furthermore, since the initial condition 
and the evolution of the mode function q p for p > 1 will be the same as the case of the flat universe inflation, we 
focus our attention on the modes with p not too greater than unity. 

To discuss the evolution in the small back reaction limit, we find it is more convenient to use the equation for the 



scalar field perturbation in the Newton gauge, Eq. (2.30), instead of working directly with (4.15). Setting 
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(4.22) 



rewriting Eq. (2.30) in terms of \ p and r]Ft., and performing analytic continuation, we obtain 

X p " + [p 2 + 2(1 - H 2 ) + d 2 Va 2 ] X p = §«0'Y ~ K^'V ■ 



(4.23) 



Also from Eq. (2.29), q p and \ p are related as 



JP — 



(4.24) 



Note that this is basically the same as Eq.( |3.11 ), except that here we are using the expansion (2.17) for q and for \, 
with the same normalization constants Afp, whereas in (3.11) we were using the Schrodinger normalization for both 
sets of eigenstates. In the limit r\ — > — oo, Eq. (4.24) reduces to 



(V -> -oo). 



(4.25) 



,±ipV 



From the above equation (4.24) and the asymptotic behavior of q p and cf>' at a — > (r) — > — oo), namely q p 
and 4>' ~ e 2?? , we find q pl ~ x p or smaller for a<l, On the other hand, the behavior of q p for a ^ 1 can be deduced 
from Eq. ( 4.15Q . Under the slow roll approximation, we find 



q p ~ e 



ld 2 V 



2 H 2 3 H 2 



(4.26) 



Hence, assuming the growing mode (the first solution in the above) dominates, we have q pl <C (cf>" /<fi')q p ~ x p . Then, 
since TL 2 -1« i/ 2 a 2 and k0' 2 = n<j) 2 a 2 , we may neglect the right hand side of Eq. ( [1.23] ) in the small back reaction 
approximation. 

Also, for the stage (Ha) 2 3> 1, we approximately have 



n p 



-K sE rX p - 



(4.27) 



Just as in the case of the flat universe inflation, this turns out to be constant in time for H 2 a 2 S> p 2 - For definitcness, 
let us explicitly show this fact and deri ve the power spectrum of 1Z C . 

The transfer function defined in Eq. (4.16) is now approximately given by 



T p « lim 

15 7)^0 



H _ 

cup 



(4.28) 



where, using Eq. ( 4.25Q , x p is th e solution of Eq. (4.23) which behaves as x p — > (2 + ip)e lv when r\ — > — oo. Now let 
us find the solution to Eq. (4.23) in the small back reaction limit. With the approximations above, and ass umin g the 
time variation of M 2 — d 2 V is small compared with the expansion time scale, the general solution to Eq. ( 4.23 ) is 



x p = aPl P (~ coth?j) + PP~, W (- cothr?), 
where — cothr/ = coshiJt, and u 1 is given by 



9 M 2 
l^Jl 2 



(4.29) 



(4.30) 



It should be noted that there are cases in which our assumptions may not be valid. For example, in a model of one- 
bubble inflation recently proposed by Linde || , the time variation of M 2 is significant at the early stage of inflation 
inside the bubble. 

The approximate solution ( [1.29 ) is valid for t < a few x H~ l , i.e., until aH is not too much greater than unity. For 
definiteness, we write this condition as aH <§C B where B is a big number, say ~ e 10 . Noticing that 



p-, lp (-cothr,) 



r(«p+i) 



(4.31) 
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we find 

X p = (2 + ip)T{ip + l)P-, ip (- cothtj) . (4.32) 

Let us consider the behavior of x p a t the stage when 1 <C H 2 a 2 <C B 2 . Note that, for p = 0(1), we hav e p 2 <C H 2 a 2 
as well. Thus we consider the stage 1 +p 2 <C if 2 a 2 <C -B 2 . At this stage, the approximate solution ( 1.32 ) is still valid 
and, using the assymptotic form of the Legendre function, the expression of x p reduces to 

(2 + ^)r(,> +1 )r(^i) e , g , 
vtt r(^' + zp + i) 

Now we consider the evolution of the background inflaton field at the stage when 1 + p 2 <C H 2 a 2 « JJ 2 . In 
accordance with our assumption, the potential V(4>) during this epoch can be approximated by 

V( C />) = V p + dV p ( C t>- C l> p ) + ^M 2 (<t>-<f> p ) 2 , (4.34) 

where 4> p is the value of <j> at a fiducial time t = t p . A convenient choice of t p is the time at which a few x (l+p 2 ) = a 2 H 2 , 
which is essentially the horizon crossing time. During this stage, one can approximate the evolution equation of the 
background field as 

4> + 3H<i> = -dV p - M 2 (0- 4> p ), (4.35) 

where H = H(t p ). The solution is 

4>-4> P = -^f + Cl e-^-^ + c 2 e-^ H ^\ (4.36) 
where c\ and ci are integration constants and 



3- ^9-4M 2 /H 2 3+^9-4M 2 /H 2 
M ^ ' A2 = 2 ' \ 4 - 67 ) 

The terms with c\ and C2 correspond to the growing and decaying solutions, respectively. Then neglecting the decaying 
solution, we can evaluate 4> during this stage as 

^ = ( dVJ e - Xl H(fi-t p ) 4 (4 38) 

H X'2 

Hence, noting that Ai = —v' + 1, we have H/ (a<j)) cx e u Ht oc x p ■ Thus TV C stays constant in time for 1 +p 2 <C H 2 a 2 <C 
B 2 . After this stage, the constancy of TV C can be shown in the exactly the same manner as in the case of the flat 
universe inflation. 

Collecting the results above, we find that the transfer function is given by 

/ ff 3 A 2 \ 2(2 + ip)r(^ + i)r(i P + i) 
s V(2a^) Al ^y t=tp V^Fr^' + ip + i) ■ 1 ' 

Notice that the factor (2aH) Xl \t=t v ~ 1 for models with M 2 « H 2 . Finally, the power spectrum of TV C at the end 
of inflation is given by 



where 5 P is a phase defined by 



H 3 \2 \ 2 2(cosh7rp + cos d p ) 



\(2aH)^dV J t=t sinh 2 7rp 



r (V + i) 



r(i^ + «p+i) 



(4.40) 



» (2-ip)T(l-ip)T(i/-ip+l)p + 



' (2 + *p)r(i + t}»)r(i/ + ip + ' 



If we set M = 0, we recover the result presented in Eq.(4.2) of hi 
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As for the supercurvature modes, although we have argued in section II] that the y are not expec ted i n thin wall 
models, they are not excluded in the general case. The normalization appearing in Eq. ( 4.11 ) and Eq. ( 4.20 ) is strongly 
model dependent. As mentioned at the end of section III, supercurvature modes can impose severe constraints on the 



parameters of the model. For a discussion of supercurvature modes in two-field models see ]8j,|9[] 

In the above discussion, we have assumed the small time variation of M as well as the small back reaction limit. 
However, as we have noted, there are models in which the former condition is violated |^J. Furthermore, it is possible 
to consider a model in which the back reaction is significant during the first few expansion times inside the nucleated 
bubble. For such a model, the spectrum at p 2 < 1 can be quite different from the one we derived above, and one should 
resort to the general expression (4.18), which typically will require numerical evaluation of the reflection coefficient and 
transfer function. Nevertheless, it is worthwhile to note that the contribution of the p <C 1 part of the scalar spectrum 
to the large angle CMB anisotropics is negligible This is due to the fact that 1/p 2 factor in the normalization 
constant at p 2 — ► is canceled by the factor p 2 coming from the harmonics \Y pim \ 2 . 



V. TENSOR-TYPE PERTURBATION 



A. Spectrum at t\r — > oo 

In our previous paper |l6[ ] , we have shown that the spectrum of the tensor-type perturbation is constrained to have 
a certain restricted form. He re we briefly review the result. The argument is parallel to the one for the scalar- type 



perturbation given in section IV A 



We expand the quantum operator w corresponding to the variable w as 

* = Y,$LKe wP (vc)Y pim + (h.c), 



(5.1) 



where w p satisfies the eigenvalue equation ( |2.39[ ). As noted in section II B , there exists no supercurvature mode in 
the tensor-type perturbation. Note also that the potential term in Eq. ( 2.39| ) goes to zero at both boundaries of the 



region C. Hence, as the two orthogonal solutions w (±y we may take those having the asymptotic behaviors as 



; , (rjc -> +oo), 



and 



(rjc -> +oo), 



As in the case of scalar-type perturbation, using the Wronskian relations we obtain 

W + \ 2 = l-\g + \ 2 7 |a_| 2 = l-|£_| 2 , 
<7_ = cr+, c+Q- + = 0. 



(5.2) 



(5.3) 



(5.4) 
(5.5) 



We also have dr/iogu^, = 2n5(p - p')<W- Thus, from Eq. ( gg7[ ), the normalization constant Af^ is determined 
as 



NT 



{p 2 + 1) 



up sinh Tip 

Finally the analytic continuation of the two independent modes to region R gives 

w r+ = e* p/2 Q + e- iprm + e-^/V^, 
W P R _ = e np/2 <7^e- %pr ' R 7 

with a constraint on coefficients, \g+\ 2 + |c-| 2 = 1. Thus we obtain 



(5.6) 



4+r 



„p i 2 



2(cosh7rp + 



(5.7) 



(5.8) 
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at rjR — > — oo. This bound is obtained without specifying details of the model. For p 2 ^> 1, the spectrum at tir — > — oo 
is model- independent. On the other hand, the spectrum for small p 2 is sensitive to the choice of a model. It should 
be noted, however, the reflection coefficient p + can be written in the form, 

p + = -\ P+ (p)\e 2 ^ (5.9) 

where |/o+(0)| = 1 and a(0) is finite. 

As mentioned in the previous section, in the case of the scalar-type perturbation, the contribution of the p 2 « 1 
modes to the CMB anisotropy is negligible. Hence we did not have to worry about the behavior of the spectrum at 
p <C 1. On the contrary, in the case of the tensor perturbation, the behavior of the low frequency spectrum does 
affect the large angle CMB anisotropy significantly. This is because there is an extra factor of 1/p 2 in the normalized 



tensor harmonics y^ ptm (see Eq. (2.43)) as compared with the scalar harmonics Y p£rn . In fact, if p + were not to 



have the behavior shown in Eq. (5.9), the CMB anisotropy due to tensor-type perturbations would have diverged. 
Hence, we must be more careful when dealing with the low frequency modes in the present case. In fact, the so-called 
wall fluctuation modes, which give an important contribution to the CMB anisotropy, are due to these modes [M. 



B. General formula for the power spectrum 

As in the case of scalar perturbations, we shall give a general formula for the power spectrum of tensor modes in 
terms of the reflection coefficient for the corresponding transfer function. 

Let us first recapitulate the evolution equation for the tensor perturbation. The equation for w p is 

d 2 



Equivalently, the equation for Q V T is 



W*-*+ +P 



a d 



w 



p _ 



0. 



dr^ a <2tr 



(5.10) 



(5.11) 



The variable w p is expressed in terms of Q V T as Eq. (2.35). Conversely, substituting Eq. (2.35) into Eq. (5.11) yields 



the expression for Q V T in terms of w p 



,2 



K 



-{aw p ). 



(5.12) 



a- a(p 2 + 1) drji 

As for scalar perturbations, we will define an auxiliary "transfer function" Tt in order to compute the power 

p I 

Tl 

k d 



spectrum for the variable Q v T ja 2 at the end of inflation . For tensor perturbations it is defined as follows, 



q-p 



lim 



iR^iend a 2 (p 2 + 1) drm 



aw 



p 



(5.13) 



where w p is the solution of Eq. ( 5.10 ) which behaves as w p — > e lpriR when t)r — > — oo. Then, taking into account the 
form of the quantum state on the lightcone inside region R, Eq. (5.7), the power spectrum of Q v T j a 2 can be expressed 
as 



( + ) |2\ 



Kl 5 



T|2 



Qt 



= 2|A/; T | 2 |7?| 2 LshTrp + KU+S 



(5.14) 



Here, we have introduced the variable U^tiL used in [p6[pi|, in terms of which the metric perturbations reads 



plm 

h l3 [QT]=Y, u( p^ 



(+) v {+)pZrn 
ij 



Thus, from Eq. ( [[!§ , we have = NJ Q p T /a 2 . 
Finally, the tensor power spectrum is given by 



( + )|2\ 



pirn I 



7?rV + i) 

np sinh ixp 



cosh np + ( Q+ -^p 



(5.15) 



(5.16) 



In the appendices we present some cases where the transfer function can be found analytically. However, for the 
remainder of this section we turn attention to the weak backreaction case. 
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C. Weak backreaction case 



Again, as in the case of the scalar-type perturbation, we omit the suffix R from all the variables and set t = tr 
throughout this subsection. We also adopt the small back reaction approximation, i.e., we assume K(j>' 2 <C H 2 . 

At the first stage when o 2 C 1, we have w p ~ e ±wn . In the small back reaction limit, this solution is valid as long 
as the inequality K<fi' 2 — K<j) 2 a 2 <C p 2 is satisfied. On the other hand, for p 2 + 1 <C H 2 a 2 , Eq. ( 5.11 ) tells us that 



I 2 a 2 , 

Qt/o. 2 becomes constant, save the decaying solution that dies off as ~ a -3 . Hence, for p 2 > 1, there exists a stage 
during which K<j) 2 <C p 2 /a 2 -C H 2 is satisfied. One can then evaluate the final amplitude of Q^/a 2 approximately at 
the t = t p when a few x (p 2 + 1) = a 2 H 2 to obtain 



Qt 



(p 2 + l 



In this case, the transfer function greatly simplifies, being given by 



(P 2 + 1) 



(5.17) 



(5.18) 



Calculating the spectrum of Q^/a 2 , we recover Eq.(4.8) of our previous paper |l6| , where the effect of the evolution 
inside the bubble was neglected. Under this assumption, the transfer function is readily given by = —kH/(j> 2 + 1). 

However, as already mentioned, the low frequency modes are very important for the tensor-type perturbation. To 
evaluate the spectrum at p 2 < 1, we perform a perturbation analysis with respect to k</>' 2 . 

First note that, once the universe is no longer exactly de Sitter, the limit 77 — ^ does not correspond to a — ► 00 
any more. To accommodate this shift in rj, we introduce a new time coordinate £ such that it satisfies <i£ = drj and 
approaches zero for a — > 00, and write the scale factor as 



-1 



a = 



Q(0 smhC 



(5.19) 



where the function Q takes care of the deviation from the exact de Sitter metric. It is assumed that Q approaches a 
constant at £ — ► — 00. Then we have 



a 
a 



coth ( 



9L 

Q 



and the Friedmann equation becomes 



Q 2 



Q' 

2sinhCcoshC— 



sinh z C ( % 



(5.20) 



(5.21) 



where the function h is defined by 



(5.22) 



The time derivative of Eq. ( 5.21 ) gives 



a 2 h 2 = 



Q 



' x O' 
-h2cothC-77. 



Q 



(5.23) 



In the limit £ — > —00, we assume that the scale factor approach that of the exact de Sitter universe. Hence we must 
have 



2e^ 2e c 

a ~^ T~ = "FT ' ' 
"0 wo 

where h = h(—oc) and Q = Q{— 00). This gives 

ho = Qoe 11 ^ =: Q e A " . 



(5.24) 



(5.25) 
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Thus the correction to the conformal time, A77, is obtained by evaluating Qq. Let us carry this out. 

To the first order in n, since Q'/Q = Of/jf), we can neglect the quadratic term in Eq. ( |5.23| ). Then it can be 
integrated to give 



9L = - 1 

Q ~ sinh 2 C 



d (> ^' 2 sinh 2 C', 



(5.26) 



where we have imposed the boundary conditi on th at sinh 2 £Q' /Q — > for £ — > 0, which means the expansion is almost 



de Sitter at late times. Inserting this to Eq. ( 5.21 ) and taking the limit £ 



2Atj _ 



= 1 + K 



d( sinh 2 CP ■ 



-00, we obtain 



(5.27) 



Hence to the linear order, we find 



An = - 
' 2 



d( sinh 2 C4>' 



12 



(5.28) 



Now let us solve for w p to the first order in n(f>' . For this purpose, we set 



e ipr >e F : 



F{n) 



(5.29) 



and insert this to Eq. (5.10). Linearizing the resulting equation, we obtain 



This can be easily solved to give 



(5.30) 



' 'rh' 2 : p 2! Ki '-v) drj' 



drj <f> e 



(5.31) 



Hence we obtain 



F = — 
2 



dr]' / dr)"(t>' 2 (ri")e 2ip ^"-T>') 



(5.32) 



where the second expression is obtained by integration by parts. 

Then the spectrum of Qt/cl 2 can be evaluated by using Eq. ( [5.17 ). Our interest here is the behavior of w p for small 
p. In particular, if the phase of w p would remain finite for p — ► 0, the tensor-type perturbation would give divergent 
contribution to the CMB anisotropy as mentioned at the end of subsection VA. Hence it is important to clarify the 
limiting behavior of w p at p — > 0. 

In the limit p <§; 1, Eq. Q5.32) reduces to 



2 



dc'(( - C')0' 2 (C) - ip I <(C - OV^C) 



(5.33) 



where we have replaced 77 with £ = 77 + Arj. The first term gives the first order correction to the amplitude of w p . 
An inspection of it shows it diverges logarithmically at £ — > 0, which implies the behavior of w p as ~ |Cl~ c where 
e = 0(n(j) 2 1 H 2 ). However, this diver gence is precisely canceled by the correction to H w ft w Q which behaves as 
- |C| e . Therefore the validity of Eq. flsHI rem ains intact. On the other hand, the second term which gives a phase 
correction which is finite at ( — > 0. From Eqs. ( [5.28] ) and (5.29), we find 



p -» |6?|e ip/3 W (C^0); /3(p) = 



dC(sinh 2 C-C 2 )0' 2 +O(rt- 



(5.34) 
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Thus the phase p(3(p) vanishes for p — > as one should have physically expected. 
Taking the above result into account, the final amplitude of Q^/a 2 is found to be 



T|2 



T|2 



cosh ixp 



(p 2 + l) 2 
p+|cos (2p{a-ff)) 



[cosh7rp — \p + \ cos (2p(a — /?))] 



2p{p 2 + 1) sinh7rp 



(5.35) 



where a and /3 are functions of p, defined in Eqs. fl5.9| ) and (5.34), respectively, and both are finite in the limit 
p —> 0. Thus the low frequency spectrum of the tensor- type perturbation is affected both by the configuration of the 
Euclidean bubble and by the evolutionary behavior of the inflaton field inside the bubble. 

We note, however, that the effect of the evolution inside the bubble will be small when the small back reaction 
approximation is valid. This comes from the fact that the integral in Eq. ( |5.34 ) can be basically evaluated over the first 



expansion times after nucleation, where sinh £ dominates and can be replaced by a, hence (3 < K(<p 2 )max/ H 2 <C 1. 
This implies the cosine starts oscillating only for p ~ 3> 1. Hence the effect of the evolution inside the bubble 
is unimportant if the small back reaction approximation is valid. On the other hand, the effe ct m ay be significant 
for a model that violates the small back reaction approximation, and then the general formula (5.16) should be used. 
This typically requires numerical solution for the reflection coefficient and transfer functions, and we have no further 
comment for such a case here. 



VI. CONCLUSION 



We presented a formalism to calculate the scalar-type and tensor-type perturbations in the one-bubble open inflation 
scenario and investigated the general features of their spectra. Our analysis is based on a single-field model of one- 
bubble open inflation, but the result for the tensor-type perturbation will be the same for a multi-field model. 

For the scalar-type perturbation, for a single-field model, we found that there is no discrete supercurvature mode 
if the wall thickness is sufficiently smaller than the wall radius. As for the continuous spectrum, we obtained the 
formula (4.4C), which reduces to the familiar expression 



<l^?| 2 ) 



9 



v 



for p 2 ^S> 1, where t p is approximately the horizon crossing time at which p 2 + 1 = H 2 a 2 . In deriving it, we adopted the 
small back reaction and assumed the small time variation of the potential curvature (i.e., the mass term). The small 
back reaction approximation is the one in which the effect of the inflaton dynamics on the cosmological expansion is 
small within a few Hubble expansion times and it corresponds to the slow roll approximation at the late stage when 
the spatial curvature of the universe can be neglected. For p 2 < 1, the resulting spectrum is largely depe ndent on 
details of a model. In particular, if the time variation of the potential curvature is large, our formula (4.40) will not 
be a good approximation. However, the effect of the very low frequency modes (p 2 <C 1) on the CMB anisotropy 
is known to be very small. Hence, we conclude that observational consequences of the scalar-type perturbations are 
rather robust except for the contribution to the low multipoles (I < a few) of the CMB anisotropy spectrum, for 
which the p 2 ~ 1 part of the scalar perturbation spectrum is important. 

For the tensor-type perturbation, it was easy to see that there is no supercurvature mode. Again, the continuous 
spectrum is found to be almost independent of details of a model as long as the small back reaction condition is 
satisfied. By taking account of the back reaction perturbatively, we found the small p behavior of the spectrum is 
affected both by the configuration of the Euclidean bubble and by the evolution of the inflaton field inside the bubble. 
For a model that satisfies the small back reaction condition, however, the latter effect is found to be unimportant. 
Nevertheless, by relaxing the small back reaction condition, it may be of interest to see if this effect gives rise to an 
additional constraint on models of one-bubble open inflation. 
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APPENDIX A: CANONICAL QUANTIZATION OF PERTURBATIONS IN ONE-BUBBLE OPEN 

INFLATION 

Here, we discuss the quantization of perturbations of a scalar field and metric in the context of one-bubble open 
inflation. The Lorentzian continuation of the Coleman-De Luccia instanton is considered as the background configu- 
ration. Following the Dirac's procedure, we reduce the degrees of freedom of a constrained system to physical ones. 
Here we follow the notation used in ]3(J . 

We start with the background metric of the form, 

ds 2 = dr 2 + a 2 (T){-d X 2 + cosh 2 x^ 2 )- (Al) 
For convenience, we introduce the unit normal vectors, 

:= (d T T = (1,0,0,0), »":=o- 1 (a K )' i = (0,o- 1 ,0,0). (A2) 
Then the metric is decomposed as 

g^u = - n^n u + a^u, (A3) 

where <7 Mt , is the metric of the 2-sphere with the radius a. Further we adopt the convention to denote the projection 
of tensors as 

V T :=V m t*, 

V n :=V^ = a- 1 V x . (A4) 

The following relation is used in the calculations. 

,. a .. tanhx ,, 

n'V = — T^n v + -a> i Vy 

q , a 

r" .„ = --n"n„ + (A5) 
a a 

The second variation of the Lagrangian for the gravitational part is given by 

4 2) = T- (- V;p'» , "' i '' + 2h^, p h"^ u - 2V^ i " + h.^) , (A6) 

OK 

where is the metric perturbation. The matter part is given by 

L%L = ^(h 2 - 2 V^) ( " + 2 
ok v a 



a „ / a \ I 



,2 



where ip is the scalar field perturbation. The explicit form of the terms in the gravitational part is given in Eq. (E.5) 
of®. 

We expand the metric and the scalar field perturbations in terms of the spherical harmonics Y — Yt m (Q) and 
consider only the even parity modes. For the metric components, we set 



l,(e) = V n f j(e)tay t(e) _ >p ir(e)tay u(e) \^ rr{e)lmy 
,v rvn / t ±J -nn J > "'nr / , ±J -nr ± ^ il tt / j j1 tt > 

u( e ) _ \ " Tr{e)£m V: .(e) _ \ " ff (e)<mv 



h (£) 
AB 



£ (w^ em Ya AB + V^ tm Y AB ) , (A8) 
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where 



Y AB = 



Y 



\\AB , 1. 



£(£+! 



(A9) 



The reality condition implies Hf m = H^~ m , where Hi = H nn , H nT , H n , H TT , H T , w, v. To keep the simplicity of 
notation, we omit the indices, (e), £ and m, unless there arises confusion. For later convenience, we list the formulas 
of the ^-integration, 



dVL YY = 1, 



dfl a AA 'Y\\ A Y\ ]A , 



dn & AA 'a BB YabYa'b* = 



ACt A-AA' ~BB' ~CC' V v 

all (7 O O 1 AB\\C* A'B'\\C 



Oil a a <J IAB\\C Y A<C'\\B> 



£(£+!), 

£(£+1) - 2 

2£(£ + 1) ' 
(£(£ + 1) - 2)(£(£ + 1) - 4) 
2l(t + 1) ' 
(£{£+ 1) -2)(£{£ +!)-&) 



U(£+l) 

It is convenient to rewrite the components having more than two of their indices projected onto the 2-sphere as 



(A10) 



h n A;B — 

h T A-B = 

hAB-n 
h-AB;T 
flAB-C 



1(1+1) TT ,l R _^X(rj W 

2a 2 cosh x a a \ a cosh \ 



£(£ + l)H n - t -^ v 



Y AB , 



w 



H T H — ( H TT — — — — 

a V o, cosh x 

1 



2a 2 cosh 2 x 
+ £(£ 
vabY 



tanhx 



oabY 



&abY 



£(£ + 1)H T - - v 
a 



Yab, 



[d x - 2tanhx]^ + ~Y AB [d x - 2tanh%]«, 
er cosh x a 
oabY 



a 2 cosh x 
w 

a 2 cosh 2 x 









'd T 2« 


w + Y A b 


'd T - 2 h - 


a 




a 



Y\\ c o A b + v Y AB \\c + -[H T - tanhx-Hn] %a0"b)c- 



(All) 



Then it is straightforward to calculate the action for the even parity modes. By using the formulas ( AlOj ) the 
fJ-integration in the action is performed to give 



d^ V~9 (4 2) + ifiL) -Jdxjdr ( L W + L«" 
Here we only demonstrate the most complicated term in the gravitational part; 

1 



(A12) 



E 



(a 2 cosh x) 3 



+8£(£ + l)a cosh 2 X [H T ~ tanh X H n ] w-2 (£(£ + l)-2)vw 



(A13) 



The matter part reduces to 



■(e) 



4«a 3 cosh x 



TT -\- TT 

nn ' tt 



a 4 cosh x 



21(1 + 1) 
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4 " I H 'r + -f^T- ( H n H r)) + 

a A cosh x J 1 



a \a I a z 



+4 



-H nT -d x X + ( H TT - 

a V 2 



)(* + /*) 



1(1 + 1) 
a 2 cosh 2 x 



H T X 



+ 



a 
i- 

a 



+ ( / + 3- ) (TrH)X + (X + fx) 1 - ^ (d x X) 



a 2 cosh x 



, +f + 3^ + f(f + 3 a 



dr 



X' 



a 
a 



(A14) 



where 



TyH — —H nn + H TT 



2w 



2 ' 



a 2 cosh x 



(A15) 



X:=%, 



(A16) 



and 



(A17) 



Here we note that there is a relation: 



b:= a - 
a 



(A18) 



which is used to remove <j) 2 from the action. As a result, the action contains k as an overall factor. We introduced / 
just for notational simplicity. In fact, it can be written in terms of b (or a), as 



t 1 d i u 



(A19) 



Hence the action is rewritten into the form that depends on the background quantity only through the scale factor a. 
This fact simplifies the calculation considerably. In addition, for computational simplicity we set 4k — 1 in the rest 
of this appendix. 

Next we define the canonical conjugate momenta by 



d{d x H\ 



{e)lr, 



(A20) 



where Pj = P nn , P nT , P n , P TT , P T , P w , P v , Px- Since the x _ derivatives of H nn , H nT and H n are not contained in the 
defining equations of the conjugate momenta, they give the primary constraint equations: 



C\ := P nn — 1(1 + l)H n — 2a cosh x sinh xH nn w + a 2 cosh 2 x 



C 2 := P nT - a 2 cosh 2 xd T (H nn + H TT )-2 



d T -2 a - 
a 



C 3 ~ P n - t(l + 1) H nn + H TT + 



a 

w = 
4 tanh \ 



d T + 2 a - 
a 



Hnr = 0, 



a 2 cosh x 



H 



0. 



(A21) 



The other components are 
2 



d x w — £(£+ l)H n — 2a cosh x sinh xH n 



2 tanh \ 



w — a 2 cosh x 



H„ 



/' =2!(l 1) ( -J) X H T -H m - 



H, 



31 



2 2 
P,„ = : 2 — [®x — 2 tanhx] w d x H TT + 2 



Pv = 



a 3 cosh x 
£(l+l) -2 1 



-<9 v v 



a 

'■ + !)- 2 



d T + 2- 



PX — uiyy/ u, won A yv x 

The Hamiltonian is defined by 



2£(£ + 1) a 3 cosh 2 x X a? cosh 2 x 
2 a cosh 2 x (c\X — aH nT ) 



» (e) = EE^ (e) M e) )- L 



(e) 



(A22) 



(A23) 



where d x H nn , d x H nT and (9 X -Zf„ are to be replaced by Ai, A2 and A3, respectively. The canonical equations of motion 
are 

H = 

P, = -i (A24) 



The consistency conditions for primary constraints are calculated as 

a 1(1 + 1) -2 



d x Ci = 2a 



%--dr+(*-- ?r , y ) --4 &a 3 COsh 2 X 

a \ a 2a 2 cosh x/ 

a a \aj a 2 2a 2 cosh x 



d T + 2/ + 3- 



-2a 3 cosh 2 x 



i? rr + 2a cosh x sinh x 



d r + 2- 



+21(1 + l)a 



a 



ff r - 2asinh 2 x#rm - 21(1 + 1) tanhx#n 



(l(l + l)-2)v 2 , . , „ 

2 h a cosh x sinh x "«j = 0, 

2a cosh x 



(A25) 



g\C 2 = 2ad T P TT + 2a 3 cosh 2 x 



9 2 + 2^ T + 2^-2 f£ 
a a \ a 



£f nT - 2Z(Z + l)a 



d T -2 a - 
a 



-2a dcosh xPw — a(P T + Px) + 4aa cosh x sinh x#nn + 4- tanhx w = 0, 



and 

d x Cz = a 



d T + 2 a - 



+1(1 + 1) 



(P T + 21(1 + l)H nT ) - 41(1 + 1) tanhx 



O t + - 

a 



2tanhx(-ffrT — H n 



a cosh x 



H n + a(P w - 2P V ) 



+ 2(1(1 + 1) -2) 



sinhx 
a 2 cosh 3 x 



(A26) 



v = 0. (A27) 



These can be solved for H nn , P v and Px- Further consistency conditions, <9 2 Ci = 0, become trivial. 

We have to mention that all the above equalities hold in the weak sense. That is, in reducing the expression we 
used the primary constraints and their consistency conditions that have been already obtained. Since the set of 
primary constraints and their consistency conditions closes, we can substitute them into the action, and remove the 
six variables, P nn , Put, Ptt, H nn , P v and Px- 

There still remain unphysical gauge degrees of freedom. In order to obtain the reduced action that contains only 
the physical degrees of freedom, we set the following gauge condition corresponding to the Newton gauge: 



2A 



b- 1 [d T + - H TT , H nT = 0, H T = 0. 



(A28) 



These gauge conditions imply the consistency conditions d x [2X + b 1 (d T + ^) H TT \ = 0, d x H nT = and 
d x H T = 0, which respectively become 



d T D x =0, A 2 = 0, D 3 



(A29) 
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where 



D 1 



D« = 



-d T - 2- + - 

a a a 



1 



1(1 + 1) 



1(1 + l)tanhx 



2a 



aP T 



21(1 + 1) 



d r -~ 



2a 2 cosh x 
(1(1 + 1) -2) ^ 

4a 2 cosh 2 x 



w + a 2 cosh 2 x 



a 



1(1 + 1) 



2a 2 cosh x 



(A30) 



From the condition d T D\ = 0, we can set D\ as an arbitrary function of \. This is due to the fact that the gauge 
condition ( A2q) does not completely fix the gauge. Therefore we must impose an additional condition. Here we choose 
Di = 0. The first two equations = and D 3 = can be solved for P T and H n . 

As the expression for D\ is rather complicated, we first consider the consistency condition for D%. The condition 
d x D 3 = determines A3 as 



A, 



u 



sinh x 
—2 tanhx ( 1 



a 

a a a 

1(1 + 1 



1 



d 2 --d T -2---- 



1(1 + 1) -4 



2 sinh 2 x 



H n — 



2a 2 cosh 2 x 
(*(* + !)- 2) 

21(1 + l)acosh 2 x 



w + a 3 coth x 



3-Or ^ 



1(1 + 1) -4 



2 sinh 2 x 



- 2 2a 2 cosh 2 x 
a 2 cothY n 
« H ^ + a(XJ a : 



(A31) 



where a(x) is an arbitrary function of x which arises due to the remaining gauge degrees of freedom. Here we simply 
choose a(x) = 0. 

Using the relations which have been already obtained, the second level consistency conditions for the first gauge 
condition, d x Di = reduces to 



a 2 + 5-d T 

a 



a 
a 



1(1 + 1) 
2a 2 cosh 2 x 



£i=0, 



(A32) 



where 



E 1 := P u 



a sinh x cosh x 



a -d T - 2 a - 
a a a 



4 Z(Z + 1) + 6 

2 



w 



-2a cothx 



3^+±-^ +i) : 6 

a a z 2a 2 cosh x 



2a 2 cosh 2 x 

(1(1 + 1) -2) 



2a 3 sinh x cosh x 



(A33) 



Again by choosing the simplest choice E\ — 0, we obtain the equation which determines P w . Furthermore d x E± = 
gives the condition, 



Ai 



sinh x cosh x 



a -d T - 2 a - 



1(1 + 1) + 6 



a 



a a^ 



-2a 2 coth x 



3^ + ± -fc^ 

i2 o„2 — „u2 



2a 2 cosh x 



ft- 



2a 2 cosh x 



1 



aa + i)-2) 
— ^ 3 — v + 

2ft 2 sinh x cosh x 2a cosh x 



(A34) 



Now we find that the set of all the constraints closes and it becomes second class. Hence we can remove the 
unphysical variables by using these constraints. Then the remaining variables are P TT , w, v and H TT . We define the 
scalar-type and tensor-type variables by 



Q 



s 

(e)lm 



a cosh x 



H, 



(e)tm 



(TF) 



= a 2 cosh 2 xHfr 



{e)tm 



(A35) 



where (TF) means the trace free part and we have used the fact that H nn is rewritten by using the set of second class 
constraints as 



2w 

Hnn SH TT + - —~2 ■ 

a 2 cosh x 



(A36) 



Then the equations for Qs and Qt reduce to 
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d x n s 

d x Q T 
d x U T 



cosh x 

-1(1 + 1) + (3 - K s ) cosh 2 x 
cosh 2 x ' 



Qs, 



-£(£ + 1)-K T cosh 2 x 



Qt, 



where II5 and IIt are defined by 

Yl( e ) em p(e)£m 

ll s .— £-T TT , 

n(e)£m T 2 

T ' := — cosh x 



£(£ + l)aH^ tm + 2 tanhx (™ (e)£m + a 2 cosh 2 xMr^ 



and the derivative operators Ks and ify are defined by 

K s := 
^ T := 



a ^ 75 — I77T~ a & 



a 75~ a 75 2 

or ar a z 



All the other variables can be written in terms of Qt/s an d ^t/s as 



H - Qs 
TT ~ 4a ' 

H QS 

tin 



4a a 2 cosh x 
2 



£(£+l)a 
H T = H nT = 0, 
4 cosh 2 x 



tanh x Qt + 



2-K T 



n, 



cosh x 

g(f + 1) + 2 
2 cosh 2 x 



Qi 



acosh v 
w = Qt H : Qs, 



tanh * n T 

cosh 2 x 



X = - 



St' 



8a6 

n T 

2 tanh x Qt § — 

cosh x 
£(£+l)-2 



1 



97 



2 cosh x 



Qi 



P -1 
ar 

P — 4 

P =^£ 



a 3 cosh 2 x 
P 2 

v e(e + i)a 3 
Px = ^-n s . 

OT 



2 tanhx Qt 



tanhx 
cosh 2 x 



cosh 2 x 



n, 



cosh x 
2 



cosh 2 x 



tanh x Qt + 2 - JsTr 



£(£ + i)-4\ n T 



2 cosh x / cosh x 



Of course, Eqs. ( A37 ) and ( A40 ) are consistent with the equations that the mode functions satisfy. 
Substituting (A4C) into the canonical form of the action 



Jdxfdr :=[d X f dr (j2J2 P t )em (dX £)lm 



h 



(e) 
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the reduced action becomes 



Jdxjdr d) = £ ( ,„ lM / +1)(£ + 2) / *C / J 



n T *r T (x T - 1) (9 x q t ) 



■-(— ^ 

2 V cosh 2 



UtKt(Kt - l)n T + Q T k T (K T - 1) + 1) + K T cosh 2 x} Qt^) 



If i 

2 Vcosh 2 



n 5 £ s n T + Q s i? s {t(t + 1) + (#s - 3) cosh 2 x } Qs 



(A42) 



APPENDIX B: SOME BOUNDS ON W IN THE THIN WALL CASE 

For the case of thin wall bubbles, the scalar field ceases to move in the true and false vacuum limits. Then the 



geometry can be approximated by a pure de Sitter one. Comparing Eqs. (3.6) and (3.9), we have 

,2 



w + w 2 



cosh 2 rj ' 



(Bl) 



where 



fj, 



H 2 



and m 2 = d 2 V/d(j) 2 (r] = ±oo). 

Solving Eq. (Bl) perturbatively, the asymptotic behavior of W is found as 



From Eq. (Bl), we can prove that 



W -> ±2 ± -n 2 e ±2v for 77 -> Too. 



W > 2 



(B2) 



(B3) 



is maintained in the false vacuum side under the condition, fi 2 > 0. The proof is as follows. If W could become 
equal to or smaller than 2, there would be a point rj = rj c at which W = 2 for the first time and W < there. But 
this contradicts with Eq. (Bl). In the true vacuum side, the mass is not necessarily large compared with H. As a 
minimum requirement, we shall assume that m 2 > in the true vacuum side. If m 2 = 0, we can solve (3.8) as 



3 + e~ 2ri 
1 + e 2 ^ 



4>' oc 2 cosh 2 rj — 3 cosh r\ sinh r\ + sinh 2 77 tanh 77 = 

cosh(r? + S) _2r, 
oc e , 

cosh 77 

where a fixed number S := > is introduced. Then Wo = W\ m 2 =0 is given by 

Wo = tanh(77 + S) — tanh 77 — 2 . 
Now, let us show W < Wo when /1 2 > —2. Namely, we prove that in true vacuum 

g = W - W < 0. 



(B4) 
(B5) 



In fact, from Eq. (B2) we have 



'-X^ 2 + Z)e- 2r i < for 77 
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Now if g = at 77 = rj c as we decrease rj from infinity, we must have g' < there. But 



9'\g=0 = 



> 0, 



cosh r^c 

which is a contradiction. Thus g < in true vacuum. Since Wq < 0, this also implies 



W 2 > Wn 2 . 



(B6) 



(B7) 



Let us denote the wall region as the region where neither (B3) nor (B7) is satisfied. Then it is straightforward to 
show that, provided that the wall is located at 77 > 0, the integral (3.17) is always smaller than tt provided that the 
wall region has a width A77 << 1. 



APPENDIX C: AN EXACTLY SOLUBLE MODEL 

Here we present a model in which both the scalar- and tensor-type perturbation equations for the evolution inside 
the bubble can be solved analytically. We follow the notation used in sections IV and V. 
Consider the case, 



2 sinh 2( 

In this case, we can solve Eq. ( |5.23 ) for Q as 



(CI) 



Q = Q (-tanhC) e , 



where 



From Eq. (5.21), we have 



2 ' '32^ 



h 2 :=^(C-f-oo) =<ft(l + 2e). 



Therefore 



In the limit £ — > —00 



Q = - 7 ij=(-tanhC) £ . (C2) 



2 > 2 
a —> —e c = —e 11 , 
Qo ho 

where in the last equality we assume that the scale factor takes the de Sitter form. Thus we find Ar] := 77 — ( as 

e 2A " = 1 + 2e. 

It is worthwhile to analyze the general behaviors of 4> and a near the nucleation point, i.e., in the limit r\ — > —00. 
If we expand 

(f> = 0o + foe 2n + 04e 4?? + • • • , 

a = A e u (1 + a 2 e 2r > + •••), (C3) 
ho 

we obtain 

a-2 = 1, 
, 3 Wo 

^ 2 = -2^' 



3G 



(> 



2- 



d 2 V 



(C4) 



where the suffix denotes a quantity at 77 — > — 00. 

The scalar field potential can be reconstructed as follows. Integrating Eq. we have 

0= -ln|tanhC|,, 
where we have set (f> — at the nucleation point. Thus 

-tanhC = e 4 ^. 

For definiteness, we assume [i < 0, which implies </> rolls toward its positive direction. Let us first examine the meaning 
of the parameter e. In terms of e, fj, is expressed as 



2e(l + e) 



We note that, from Eqs. (C3) and (C4), we have 



3dV , 
^^ 1 + 2 ^ 



which implies 



dVo 
V 



/32Ke(l + e) 
9(1 + 2e) 2 



where Vo and Vq' are the values of the potential and its (^-derivative, respectively, at the nucleation point. Furthermore, 
expanding cf> to the next order in e 2 ^, we find (j>4 = in the present case. Hence 

A/ 2 := d 2 V = 2hl . 

Thus the curvature of the potential is relatively l arge at th e nuc leation point. 

Now we reconstruct the potential. From Eqs. ( [5.22| ) and (5.21) together with Eq. (C2), the potential V is given by 



V = -Q 2 

K 



(3 + 2e)(l + e) + e(l-2e) tanli2c 



Vo 



3(1 + 2e) 
3(1 + 2e) 



exp 



2ne 



1 + e 
(3 + 2e)(l + e)exp 



(3 + 2e)(l + e) + e(l-2e)exp 
- e(l — 2e) exp 



2 k 



2k6 



1 + e 



e(l + er 
2k(1 + e) 



(C5) 



Assuming e<l, the potential becomes very flat for > \fejn. 



Now let us consider the scalar-type perturbation. The equation (4.15) for q p in this model becomes 

d 2 AA 2 



d( 2 sinh 2 2C 



Setting z = e 4< >, this equation reduces to 



d 2 1 d 



A 2 



P 



A 2 



q p = 0. 



dz 2 z dz (z — l) 2 2(2—1) 16z' j 
This can be solved in terms of a hypergeometric function. The solution is 



q p = 0. 



^ = e *P(C+ATj)( 1 _ e 4C)l+e 2 .F 1 



l + e,l + e + | ;1 + | ;e 4 < 



(C6) 
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where we have already choosen the boundary condition q p 
background and of q p as £ — > 0, we find 



as 77 



ftO ^ 

4 s 



-00. Computing the behaviour of the 



(C7) 



JpAv r(l + 2e)r(l + »p/2) 

r(i + e)r(i + e + »p/2V w 



XI + e)r(l + e + ip/2) 
from where the scalar transfer function is found to be given by 



n-v 



hoe ipAv 2 -2e-l/2 x / K 



(e + l)(2e+l) r(l + 2e)r(l + ip/2) 



r(l + e)r(l + e + ip/2) 



(C8) 



For the tensor- type perturbation, the equation (5.10) for w turns out to be exactly the same as the one for q in 
this model. Therefore one has w = q, and thus 



{fip _> e ipv — pip(C+&v) 



(C--00), 



ipAv T(l + 2e)r(l -Hp/2) 

r(l + e)T{l + e + ip/2) [ ^ 



(C-0). 



(C9) 



In particular, the phase (3 introduced in Eq. (5.34) of section |V C| is found as 



(3 — An + - arg 
P 



r(i + *p/2) 



r(i 



ip/2) 



(l--) £ (e«l) 



(CIO) 



The last line agrees with the analysis in section V C under the small back reaction approximation. Finally, the tensor 
transfer function is given by 



T (p 2 + l)T(l + e)T(l + e + i P /2) 



(Cll) 



APPENDIX D: EXACTLY SOLUBLE EXPONENTIAL HT MODEL 

In ]3l| ] an analytically soluble exponential model in the context of the Hawking- Turok Open inflation was considered. 
In this model the potential is given byQ 

v= '¥Ii e V^t (Di) 

The exact analytical solutions for the background are [|T) 

y/2 tanh 1/2 r) C 

a (Vc) = 77 r ( D2 ) 

Hq cosh r/c 

Hvc) = -\ 7T" IntanhTyc, (D3) 
V 2k 

where rjc — — ln(tanr/2). The range of r\c is (0,oo). In this model the reflection coefficient g + can be analytically 
computed for both scalar and tensor perturbations, so a closed expression for the power spectrum can be given (notice 
that the transmission coefficient cr_ vanishes because the singularity acts as a reflecting boundary). 

We should mention that perturbations in a Hawking- Turok model with linear potential has been considered in p3[. 



5 In fact, this potential is a particular case of the potential discussed in appendix C with e = 1/2. There only the evolution of 
the modes after nucleation is considered, whereas here we consider the normalization of the primordial spectrum. 
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1. Scalar power spectrum 



In order to compute the power spectrum of the scalar perturbation on the comoving hypersurface 1Z C , we follow the 
formalism described in section IV. We need to find the reflection coefficient g + and the special solution q p for this 
model. 

To com pute g + in Eq. (4.4), we need to solve the equation for the perturbations outside the light-cone. The 
equation ( E-jj ) for q turns out to be 



pi i 



sinh 2rjc 



q p =p 2 q p . 



(D4) 



The general solution of this equation can be expressed in terms of a hypergeometric function (just taking Eq. ( JC6| ) 
with e = 1/2 and C = ~Vc): 



3 3 + ip ip 4no 

^ ? r-. ? I" n 1 C 



(D5) 



and its complex conjugate. This particular solution tends to e~ tpric as rjc — > oo and diverges as t)q 2 when 77c* ~~ * + . 
Then, the normalized solution which corresponds to a wave coming from rjc — 00 which bounces against the singularity 
at 77c — * + is 



iq 1 . 



r(i + i P /2)r(3/2 - ip/2) . 



r(l -zp/2)r(3/2 + zp/2) 
which is everywhere regular. The reflection coefficient g+ can be read from 

r(l + ip/2)r(3/2 - ip/2) 



r(i - ip/2)r(3/2 + ip/2) ' 



(D6) 



(D7) 



Notice that in this case a complete set of modes is given by q p + with p > 0. 

The modes are continued inside the light-cone by means of the analytical continuatio n rj c = ~r\B. — mt/2. Their 
evolution inside the bubble is just then given by the analytical continuation of Eq. ( |D6| ). In fact, the analytic 
continuation of q p is just the solution q p we need to compute the transfer function 



j P = e < P"*(l-e 4 '' H ) 3 / 2 2 Fi 



3 3 + t P -1 | W. c i nR 



2 2 

Now we compute the assymptotic form of q p and of the background as an — * 00 . As t]r goes to 0, we find 

- P , 1 r(i + » P /2) 

^ ^(-%) 1/2 r(3/2 + zp/2)' 

IT 1 

2kt?_r ' 
V2 1 
^o(-w) 3 / 2 ' 



Using the results above, the scalar transfer function (4.16) for this model turns out to be 

':-!/,- F(l + ip/2) 



T$ = H 



ir (p 2 + l)r(3/2 + zp/2)' 



Using Eq. ( 4.18Q , we finally find the primordial power spectrum for the scalar curvature perturbation, 

1 



(\K\1 



3kH 2 



7T (p 2 +4)(p 2 + 1)' 



(D8) 

(D9) 
(D10) 
(Dll) 

(D12) 

(D13) 



The power spectrum of temperature anisotropies predicted for this model due to scalar perturbations is shown in 
Fig. @. 
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FIG. 6. The spectrum of temperature anisotropics predicted for the exponential Hawking- Turok model described in ap- 
pendix^. We have choosen the cosmological parameters h = 0.70, f2 m = 0.35, £7b=0.05, Q\ — 0, Yhz = 0.24, and 7V„ = 3.04. 
We show the individual contributions from the scalar (S) and tensor (T) modes, as well as the total contribution. Since inflation 
never ends for an exponential potential, we have artificially truncated the potential so that the number of e-foldings corresponds 
to a present value of the density parameter fio = 0.4. Notice that the spectrum is finite in spite of the singular nature of the 
background. 



2. Tensor power spectrum 



For the present model, the equation ( [2.39 ) for the tensor perturbation variable w p coincides with the one for q 



w ' 



sinh 2rjc 



w p = p io f 



(D14) 



The reflection coefficient in Eq. (5.2) for the tensor perturbation is then the one computed in the previos section, Eq. 
(D7), and the solution w p coincides with q p , Eq. (D8). 



Collecting these results, the tensor transfer function (5.16) for this model is given by 

2kH 



n-v 



r(l + ip/2) 
/2lF (p 2 + l)r(3/2 + ip/2) 



and the power spectrum by 



(\u, 



(+) |2\ 









Qt 




a 2 



8KiJ 2 
7T(p 2 + l) 2 



(D15) 



(D16) 



The power spectrum of temperature anisotropies predicted for this model due to tensor perturbations is shown in 
Fig. |. 
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